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SVDS PLUME IMPINGEMENT MODELING DEVELOPMENT 


SENSITIVITY ANALYSIS SUPPORTING LEVEL B REQUIREMENTS 


By P. B. Chiu, D. J. Pearson, P. M. Muhin*, and P. B. Schoonmaker, 
McDonnell Douglas Technical Services Co., Houston, Tx.; and R. J. Radar, 

MDAC, Huntington Beach, Calif. 


1. SUMMARY 

This report describes a series of sensitivity analyses (trade studies) 
performed to select features and capabilities to be implemented in 
the new Plume Impingement Model being developed under Task Order 
D0717 of Contract NAS9-14960. 

The objective of these studies was to ensure that this model will 
provide the performance required by the user community, but not be 
"loaded down" with features and capabilities for which there is no 
real need. A subordinate objective was to demonstrate that the cor- 
rect decisions were made in defining the model capabilities, and 
thus build confidence in the model among the user community. The 
recommendations of this report are to be used to establish a final 
set of Level B Requirements for the Plume Impingement Model. 

Sensitivity analyses were performed in study areas pertaining to 
geometry, flowfield, impingement, and dynamical effects. Recom- 
mendations based on these analyses are summarized in Table 1-1, 
reproduced from Section 5 for convenience. 
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TABLE 1-1 SUMMARY OF MODELLING RECOMMENDATIONS 


MODELLING AREA 

RECOMMENDED IMPLEMENTATIONS 

4.1 

Orbiter/Jet Geometry 

• Include shadowing caused by the Orbiter mol'dline, cabin, OtlS-pod bulkhead, and tall. 

• Do not consider reflections of the plume. 

4.2 

Object Geometry 

• Include some "first order" simplification In the object modelling. 

• The mesh si^e should be variable, and on the order of 10x10 In the near field. 

• Include back-side shadowing, and shadowing of one payload element by another payload subshape. 

4.3 

Flowfleld 

l 

• Model the flowfleld as patched source flow with radial streamlines. 

• Neglect effects on the flowfleld due to nozzle scarfing and ambient atmosphere effects. 

• Include a "practical boundary" test. 

• If simultaneous firings of adjacent engines occur frequently, implement an ."equivalent engine” 

formulation. . • ■ 

4.4 

Impingement 

• Calculate Impingement pressure in the continuum regime using a modified newtonian formulation. 

• Calculate Impingement pressure in the free-molecule regime using thermal .accommodation analysis; 
use approximations for large impingement angles. 

. * * O 

• In the transitional regime, use a linear weighting function, based on the sin function, to 
combine continuum and free-molecule results. 

1 

4.5 

Dynamical 

Effects 

• Neglect plume Impingement on the RMS for integrated simulations. Use off-line "stand-alone" 
analysis to Identify gross RMS operating constraints. 

4.6 

Contamination 

To be supplied. 

I 


“Additional NASA Input required for final decision. 




2. INTRODUCTION 


It has been clear for some time now that a plume impingement modelling 
capability is required to conduct meaningful studies of payload proxi- 
mity operations and payload handling operations. This is because the 
effects of RCS thruster plume impingement— 'both payload disturbances 
and contamination effects (see Figure 2-1)— are potentially the domi- 
nant factors driving the development of techniques for these operations. 


It also appears likely that new modelling software will have to be de- 
veloped to support the requisite mission analysis and mission planning 
studies. Existing plume-modelling software falls into two general 
classes: (1) hi gh- fidelity programs which are too large, too slow and 

too unwieldy to be integrated with total on-orbit simulation programs, 
and (2) simple models which may not provide sufficient performance to 
satisfy the requirements of the anticipated Mission Planning and 
Analysis Division (MPAD) users. What is desired is a model providing 
high accuracy, fast execution, and small storage requirements. 
Reference 1 documents the final Level A requirements and structure 
for such a model, which is being developed under Task Order D0717 
of contract NAS9-14960. 


The objectives of this Plume Impingement Model Sensitivity Analysis 
Study-Report are to establish the tradeoff factors between "cost" 
(in terms of storage, execution time, complexity) and performance 
(accuracy, flexibility), and thereby support level B decisions on 
modelling approaches and features. 
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FIGURE 2-1 PLUME IMPINGEMENT PROBLEM 


3. STUDY METHODOLOGY 


This section describes the sensitivity analysis study breakdown, and 
the analysis tools and basic data which were used by the MDTSCO/MDAC- 
West study team in performing the individual studies. 

3.1 Study Breakdown 

The basic organization of the Sensitivity Analysis Study task follows the 
Level A function/module breakdown of the Plume Impingement Model (Reference 
1). Within each study area (e.g., Orbiter geometry, flowfields), model- 
ling issues or alternate approaches are identified. No modelling issues 
are identified in the contamination area, since implementation of this 
module is not yet defined in detail. A dynamics study area was also esta- 
blished, even though vehicle dynamics are not part of the plume model 
itself, to explore implications of vehicle dynamics upon plume model require- 
ments. Figure 3-1 shows the sensitivity analysis study breakdown super- 
imposed upon the Level A function/module breakdown. 

For each study area, individual studies were conducted to resolve the issues 
or decide among the alternate approaches; these individual studies are de- 
scribed in ‘Section 4. Finally, the results of these independent studies are 
correlated, and potential interactions are considered. The results are 
integrated and presented in Section 5. 

3.2 Analysis Tools and Basic Data 

In the sensitivity analysis study, certain data and programs were used 
repeatedly. Data concerning the RCS thrusters and, scientific payloads 
are listed here, along with general information pertaining to flowfleld 
impingement, and contamination programs. 

3.2.1 Orbiter RCS Data 


The RCS thrusters (Reference 8} are Marquardt #R40 engines generating 870 
pounds of thrust with an MMH/N^O^ propellant system. These engines can 
either fire continuously or operate in a 40 msec pulse mode. The oxidizer/ 
fuel ratio is 1.6 and the chamber pressure is 152 psi. 
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FIGURE 3-1 PLUME MODEL LEVEL A RJHCTIOH/M03UIE BREAKDOWN (FROM REFERENCE 1), SHOWING SENSITIVITY. ANALYSIS STUDY AREAS 










The thrusters have bell nozzles with various amounts of scarfing (nozzle 
extension) determined by the Orbiter mold line; see Section 4.3.2 . For 
the basic (unscarfed) engine, the assumed lip angle is 16-0 degrees, the 
exit radius is 4.8 inchds, and the expansion ratio (exit area: throat- 
area) is 22:1. In this report, unless otherwise stated, the abbve basic- 
engine parameters are used. 

These RCS engines are located in forward and aft modules on the Shuttle 
Orbiter, as shown in Figure 3-2 (Reference 8). The thruster identifica- 
tion number, thrust direction, and assumed location in station coordinates 
are presented in Table 3-1. 

3.2.2 Payload Data 

The two payloads of greatest interest for this study are the Long Duration 
Exposure Facility (LDEF) and the Multi-mission Modular Spacecraft-Solar Max- 
imum Mission (MMS/SMM). 

LDEF: The NASA/Langley LDEF is 30 feet in length; its actual cross-section 

is that of a twelve-sided regular polygon having a major diameter of approx- 
imately 14 feet (Figure 3-3). The primary structure consists of aluminum 
ring frames and longerons. This skeleton weighs approximately 7,000 pounds. 
Seventy-two experiment trays are mounted in the bays formed by the rings 
and longerons, adding about 7,000 pounds more in weight. LDEF is gravity- 
gradient stabilized and includes an 8-inch spherical viscous magnetic dam- 
per with a damping coefficient of up to 2 ft-1bs/(rad/sec) (Reference 2). 

•Unless stated otherwise in this report, LDEF is assumed to be a cylinder 

(30* x 14’ D), weighing 15,000 pound's, with no damper. The mass center is 

assumed to be at the geometric center. The minimum principal' moment of 

inertia is assumed to be 20,000 slug-ft 2 with the remaining two moments of 

2 

Inertia equal to 50,000 slug-ft . 

j^S-SMM : The first of the Goddard Space Flight Center multimission modular 

spacecraft is the Solar Maximum Mission configuration, shown in Figure 3-4. 
The service section of the payload is the basic MMS structure. This roughly 
triangular section contains replaceable subsystem modules, including an 
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FIGURE 3-2 RCS THRUSTER NUMBER, DIRECTION OF THRUST, AND LOCATION 






TABLE 3-1 RCS THRUSTER IDENTIFICATION AND LOCATION 


Forward RCS Thrusters 

Thrust 
Direction ■ 

Thruster 

Identification Number 

Thruster Locations in Station Coordinates* 1 , 
in. 

X 0 

Y 0 

Z 0 

Old 

Hew 

-X 

■RMl 

F2F 

327.277 

14.654 

392.955 



F3F 

327.383 

0.0 

394.450 



F1F 

327.277 

-14.654 

392.955 

+Y 


F1L 

362-. 671 

-55.631 

373.728 



F3L 

364.708 

-58.579 

359.250 

-Y 


F2R 

362.671 

55.631 

373.728 


144 

F4R 

364.708 

58.579 

359.250 

+Z 

135 

F2U 

350.925 

14.394 

399.588 


125 

F3U 

350.917 

0.0 

400.818 


115 

F1U 

350.925 

-14.394 

399.588 

-Z 

136 

F2D 

333.840 

49.814 

372.350 3 


116 

FID 

333.840 

-49.814 

372.350 a 


146 

F4D . 

348.440 

54.839 

373.566 a 


126 

F3D 

348.440 

-54.839 

373.566 a 

-Z (Vernier) 

158 

F5R 

324.350 

53.830 

357.900 3 


157 

F5L 

324.350 

-53.830 

357.900 a 



Aft RCS 

Thrusters 



+X 

' 331 

R3A 

1555.293 

137.00 

473.058!' 


311 

R1A 

1555.293 

124.00 

473.058:* 


231 

L3A 

1555.293 

-137.00 

473.058:* 


211 

L1A 

1555.293 

-124.00 

473.058® 

+Y 

243 

L4L 

1516.00 

-135.21 

459.00 


223 

L2L 

1529.00 

-135.21 

459.00 


233 

L3L 

1 1542.00 

-135.21 

459.00 


213 

UL 

1555.00 

-135.21 

459.00 

-Y ' 

344 

R4R 

1516.00 

135.21 

459.00 


324 

S2R 

1529.00 

135.21 

459. 00 


334 

R3R 

1542.00 

135.21 

' 459.00 


314 

R1R 

1555.00 

135.21 

459.00 

+Z 

245 

L4U 

1516.00 

-132.00 

iHH 


225 

L2U 

1529.00 

-132.00 

480.50 


215 

111) 

1542.00 

• -132.00 

480.50 


345 

R4U 

1516.00 

132.00 

480.50 


325 

R2U 

1529.00 

132.00 

480.50 


315 

R1U 

1542.00 

132.00 

480.50 

-z 

246 

L4D 

1516.00 

-111.945 

437.403 


226 

12D 

1529.00 

- 111.00 

f wTiliMraB 1 


236 

L3D 

1542.00 

-110.055 

442.597 


346 

R4D 

1516.00 

111.945 

437.403 


326 

R2D 

1529.00 

111.00 

440. 00 c 


336 

R3D 

1542.00 

110.055 

442.597 C 

-Z (Vernier) 

257 

LSD 

1565.00 

-118.00 

455.44 


358 

R5D 

1565.00 

118.00 

455.44 

-Y (Vernier) 

357 

R5B. 

1565.00 

143.38 

459.00 

+Y (Vernier) 

258 

L5L 

1565.00 

-143.38 

459.00 


?Csnted 37° outboard in the Y— Z plane. 
"Canted up 10° in the X-Z plane. 

C Canted aft 12° in the X-Z plane and 20° 


d Thr 


uater mount attach points. 


outboard in the Y-Z plane. 
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FIGURE 3-3 LDEF PAYLOAD CONFIGURATION AND NOMINAL ON-ORBIT ORIENTATION 
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attitude control system of reaction wheels and magnetic torquers. The max- 
imum capability of the momentum wheels to absorb external disturbance torques 
and maintain control of the spacecraft is approximately 12 Ib-ft-sec 
(Reference 3). 

Unless stated otherwise in this report, the MMS/SMM payload is modelled as 
a cylinder of 14' length and 7.2' diameter with two solar panels 8.7 1 on 
a side (see Figure 3-5). Its mass is assumed to be 4000 pounds, with prin- 
cipal inertias being 2000, 2000, and 680 slug-ft 2 . The center of mass is 
assumed 10" from the geometric center of the cylindrical body. 

3.2.3 Flowfield Programs 

Two computer codes are used in sequence to provide plume flowfield data re- 
quired for the sensitivity analysis study: (1) the one-dimensional equilibrium. 
(ODE) program (Reference 4) developed at NASA/Lewis, and (2) the method of 
characteristics (HOC) program developed by Lockheed/Buntsville (Reference 5). 
Both programs have been developed, improved, and extensively used for more 
than ten years. They are considered the "standard" programs in plume char- 
acterization by the industry. 

ODE : Basically, this program calculates frozen or equilibrium chemical com- 

positions, species concentrations and thermodynamic properties for a com- 
plex chemical system in a variety of gasdynamic flow problems under the 
assumption of a quasi -one-dimensional flow. 

The rocket performance calculation option of ODE was used to obtain the 
thermochemical properties required by the MOC program. To obtain these data 
for a rocket nozzle, ODE requires information on chamber pressure, pro- 
pellants, and oxidant/fuel (O/F) mixture ratio (see Section 3.2.1). In 
addition, either a prescribed pressure drop or area ratio is needed in 
the rocket performance calculation. 

HOC : This program was developed for calculating a steady, two-dimensional, 

supersonic, and compressible flow in a rocket nozzle and/or plume Including 
the effects of real gas and 0/F gradients. 
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To run this program, the thermodynamic properties of a real gas mixture, 
and a description of nozzle geometry and external pressure boundary con- 
ditions must be input. The program then generates the flowfield and ther- 
modynamic properties along "left-running characteristic" lines inside the 
nozzle and plume (see Figure 3-6). To provide the flowfield data for the 
plume impingement (PUMP) program (to be discussed in the following section), 
we have used the nozzle exit plane start option of the MOC program with the 
real gas properties obtained from ODE* and nozzle parameters as listed in 
Section 3.2.1. The nozzle flow was assumed to expand into a vacuum. 

No shock waves were simulated in the flowfield, due to the exit plane start 
option; however, both inviscid and viscous flowfield data (i.e., including 
boundary layer) were generated by the program. These flowfields were cal- 
culated out to 250 feet in the axial direction and 200 feet in the radial 
di recti on. 

Although the MOC program is very versatile code, it becomes inefficient 
and less accurate as it tries to compute beyond the continuum regime of 
the flowfield. (See Section 4.3 for methods of providing flowfield data 
at greater distances.) 

3.2.4 Impingement Programs 

Two computer programs used in the sensitivity analysis study for computing 
plume impingement effects are (1) the PLIMP program developed by Lockheed/ 
Huntsville (Reference 6), and (2) the Quick-Look Disturbance Program 
(QLDP) developed by MDTSCO (Reference 7). Like ODE and MOC, PLIMP has 
also been developed, improved and extensively used for more than ten years. 

The QLDP program was used mainly for a quick-look "order of magnitude" analysis. 

PLIMP : This program computes forces, moments, and heating rates on sur- 
faces immersed in or subjected to a one-rocket exhaust plume environment. 
Program capabilities include the treatment of continuum, transitional, and 
free molecular flow regimes for arbitrary body shapes. In addition, body 
shadowing effects can also be calculated by the program. 




FIGURE 3-6 GENERATION OF FLOWFIELD PROPERTIES ALONG 

LEFT-RUNNING CHARACTERISTIC LINES USING MOC 


15 


To run this program, plume flowfield data generated with the MOC program, 
as well as descriptions of the engine, reference, and subshape coordinate 
systems and subshape geometries, must be input. The PUMP program first 
reads in the MOC-generated raw flowfield data from a tape, then rearranges, 
orders, and stores these data for more convenient access. Impingement 
calculations are then performed on each elemental area resulting from sub- 
division of the subshape (see Section 4.2.2). The resulting forces and 
moments are output in a specified coordinate system. 

PLIMP requires fairly long computer time per run case because of its com- 
plexity and versatility. 

QLDP ; QLDP is a quick-look disturbance program developed by MDTSCO for 
determining "order of magnitude" plume impingement disturbances. It is 
used for offline analysis of approach trajectories, computing plume-in- 
duced payload dynamics. 

The payload and Orbiter are modelled in only two dimensions; a plane of 
symmetry containing all jets is assumed. The Orbiter- relative payload 
initial state is input, along with time histories of payload state and 
thruster firings. Impingement forces and moments are computed using a 
MDTSCO-developed curve fit to MDAC-West data (Section 4.3.1) and the new- 
tonian impingement formulation (Section 4.4), neglecting shadowing 
(Sections 4.1.T, 4.2.3). State perturbations due to these forces and 
moments are computed and propagated for each time interval. All motion 
is three degree-of-freedom: two linear coordinates and one angle. 

3.2.5 Contamination Programs 

(This section is to be supplied at a later date.) 
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4. , INDIVIDUAL SENSITIVITY ANALYSIS STUDIES 


Individual sensitivity analysis studies are independently documented in 

the following subsections. Each individual study deals with a specific 

modelling issue in a specific area of the Level A Plume Node! Function/ 

Module Hierarchy (Section 3.1). 

* , 

A standardized structure is used for the documentation of each individual 
study, in each 4. subsection: 

The question - What is the question to be answered or what are 
the modelling alternatives to be considered? Some background 
information will often be included, relating to physical phenomena, 
operational aspects, or other considerations. 

4._._.2 Tradeoff factors /considerations - How will the alternatives be 
implemented? What will be the qualitative "cost" in terms of 
complexity, storage, execution time, etc.? How will the alterna- 
tives affect the performance of the model (i.e., accuracy,, 
flexibility)? 

4._.„.3 Analysis approach - This section explains the theory behind the 
alternatives and, if the problem is to be solved analytically, 
how the solution will be obtained. For those problems using 
computer programs the type and number of runs made is included 
here along with any assumptions made. If any payloads are to 
be used in the analysis, their position and orientation are also 
given. 

4._._.4 Analysis results - The results of the mathematical and/or com- 
putational analysis are presented and evaluated in terns of 
cost and performance. 

4._._.5 Preliminary conclusions and recommendations - Based on evaluation 
of the analysis results, conclusions and recommendations are 
made. These are preliminary since they depend only on the indivi- 
dual study under consideration; they may be revised when the in- 
teraction between the individual studies is considered in Section 5. 
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4.1 Orbiter/RCS Geometry Modelling Approaches 

The studies documented in this section deal with alternate formulations 
of the Orbiter/jet geometry module, ORBJET. The following functions are 
allocated to ORBJET (Reference 1): 

• Determine the jet classification (primary, vernier, scarfed,. . .), 
using the jet tag and a table look-up. 

• .Determine the coordinate transformation from the object-fixed 
axes to the thruster axes for the proper jet. 

• Determine, if necessary, Orbiter shadowing parameters (size and 
location of Orbiter surfaces which may shadow parts of the 
payl oad ) . 

• Determine, if necessary, Orbiter reflection parameters (size and 
location of Orbiter surfaces which may reflect plume exhaust onto a 
payl oad). 

The following subsections describe studies which were conducted to answer 
the following questions: 

• Is the capabi 1 i ty- to simulate Orbiter shadowing of RCS plumes 
required in the model? 

• Is the capability to simulate Orbiter reflection of RCS jplumes 
required in the model? 

Basic Orbiter and RCS data used in the sensitivity studies are defined 
in Section 3.2.1. 
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4.1.1 Orbiter Shadowing of RCS Plumes 

4. 1.1.1 The Question 

Is it necessary to include an Orbiter shadowing capability in the plume 
impingement model? 

Two types of Orbiter shadowing of RCS plumes may be considered. One 
of these is due to the moldline (Orbiter skin contour) around a 
thruster port and the other is due to the intruding Orbiter structural 
shapes,, such as the cabin, 'wings, tail and bulkheads. Both types of 
Orbiter shadowing can block the "line-of-sight" from a point on a 
payload or RMS arm to a thruster (similar to the payload subshape/sub- 
shape shadowing considered in Section 4.2.3). 

4. 1.1. 2 Tradeoff Factors/Considerations 

Implementation : To model Orbiter shadowing, we must first identify those 
parts of the Orbiter structure in the vicinity of RCS thrusters which 
may cause shadowing of a plume during payload proximity operations, 
then generate parametric data representing their geometry. 

To implement moldline' shadowing, the model would represent the Orbiter 
skin by the tangent plane at the thruster port, and implement the test 
with a simple vector dot-product operation. 

Intruding objects would be represented as either circular or polygonal 
flat plates, specified in Orbiter station coordinates. For a particular 
thruster and object, the model would then generate the thruster-referenced 
parameters of the "shadow zone" created by the object: a cone for a circular 

plate, a "pyramid" for a polygonal plate. Each element is then tested to 
see if it falls inside any of the shadow zones for that thruster. 
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"Cost" Factors : Implementing Orbiter shadowing would require extra 

storage and increase complexity in programming. Additional execution time 
is also required for the shadowing tests themselves. However, these tests 
may save some execution time, due to bypassing impingement computations — 
particularly for close-in operations. 

Performance Factors : Neglecting the Orbiter shadowing effects would be 

expected to cause significant force/moment errors for payloads at 
close ranges. 

4. 1.1. 3 Analysis Approach 

Three different Orbiter substructures were considered for assessing Orbiter 
shadowing effects (Figure 4. 1.1-1): (1) the Orbiter cabin for the forward 

+Z jet shadowing, (2) the Orbiter OMS-pod bulkhead for the aft +Z jet shadowing, 
and (3) the Orbiter moldline for the forward +Y jet shadowing. 

For modelling Orbiter shadowing, both the cabin and the OMS-pod bulkhead 
were represented by circular plates with the plate diameters taken equivalent 
to the respective substructures in shadowing. The moldline was represented 
by a flat plate placed parallel to the nozzle exit plane. Figure 4. 1.1 -2 
depicts the various Orbiter shadowing representations. 

Various objects including LDEF, MMS/SMM, and RMS were included in the 
Orbiter shadowing study. LDEF, MMS/SMM, and RMS were realistically placed 
at various locations and orientations for the cabin- shadowing effect, 
while only LDEF was considered for the moldline and the bulkhead shadowing 
effects. The individual object location and orientation with respect to 
each RCS engine are also shown in Figure 4. 1.1-2. 
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FIGURE 4.1. 1-1 


jets and orbiter substructures analyzed in the 

ORBITER SHADOWING STUDY 
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(a) Cabin Shadowing of Fwd +L Jet on LDEF, MMS, and RMS 



(b) Moldline Shadowing of Fwd +Y Jet on LDEF 


y 




(c) Bulkhead Shadowing of Aft +Z Jet on LDEF 


FIGURE 4.1. 1-2 SCHEMATIC DRAWINGS OF VARIOUS ORBITER SHADOWING 

REPRESENTATIONS AND OBJECT LOCATIONS AND ORIENTATIONS 
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Raking use of the PL.IMP program (Section 3.2.4), impingement calculations 
were performed for all the Orbiter shadowing combinations discussed in 
the previous paragraph. Each run was made with and without the Orbiter 
shadowing option enabled. Both inviscid and viscous flowfields 
were provided by MOC (Section 3.2.3). The inviscid flowfield was used in 
the cabin and bulkhead shadowing cases, while the high angle viscous flow- 
field was used in the moldline shadowing case. 

Modified Newtonian impingement pressure (Section 4.4:1) was used in all 
the runs made; the mesh fineness used was between 10 x 10 and 20 x 20 
(Section 4.2.2) for all the object subshapes. 

4. 1.1.4 Analysis Results 

Table 4. 1.1-1 provides case-by-case performance data for the runs made 
to analyze Orbiter shadowing. Note the substantial errors occurring for 
some test cases -- over 100% in force magnitude and 10° in force direction; 
up to 50% inmoment magnitude and 5° in moment direction. 

Typical computer time in terms of Uni vac 1108 Central Arithmetic Unit 
was roughly 25 seconds per run. Furthermore, all the runs with shadowing 
took slightly less time than those without shadowing. 

4.-1.1.5 Preliminary Conclusions and Recommendations 

Based upon the above study, we have found that neglecting Orbiter shadowing 
effects can cause significant errors in both force and moment calculations. 
Furthermore, inclusion of the shadowing effects may actually save computer 
time. 


23 




TABLE <.1.1-1 PERFORMANCE VARIABLES FOR ORBITER SUBSHAPE SHAOCMMG 
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* L denotes LBEF 
H denotes KJ5/SIW 
R denotes RMS 

h See Figure 4.1. 1-2 

c With respect to the results with Orblter shadowing 




































Therefore, we recommend that the Orbiter shadowing capability be imple- 
mented in the plume impingement model. 
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4.1.2 Qrbiter Reflection of RCS Plumes 

4. 1.2.1 The Question 

Is it necessary to account for the effect on the payload due to the reflection 
of the RCS plumes from Orbiter surfaces? Since the reflection mechanism is 
different depending on the impingement regime (free molecular, transitional, 
or continuum), the question must be answered for each regime. These impingement 
regimes are discussed in Section 4.4. One Orbiter surface which could 
reflect an RCS plume enough to effect a payload is the tail, which potentially 
could reflect an aft +Z RCS engine plume. 

4. 1.2. 2 Tradeoff Factors/Considerations 

To assess cost and performance factors, consideration of how plume reflection 
would be implemented in the model is necessary. 

Impl ementation : To Implement plume reflection In the model, a technical effort 
would be required to analyze and develop techniques to calculate the plume 
reflection and the effect of this reflection on the payload. The developed 
techniques and associated equations and logic would then be incorporated into 
the model . 

"Cost 11 Factors : The costs of implementing plume reflection are: (1) the cost 

associated with developing the necessary techniques and implementing the tech- 
niques into the model, and (2) the associated storage and execution time. 

Performance Factors : Since the effect on the payload due to plume reflection 

from the tail is expected to be minor, the increase in accuracy in the calculation 
of the plume impingement effects including the effects of plume reflection 
should be minimal . 

4 . 1 . 2 . 3 Analysis Approach 

The impingement Knudsen number (ratio of the molecular mean-free-path in the 
plume, A m , to the characteristic dimension of the impinged surface, L s ) was 
calculated for the case of an aft +Z RCS plume impinging the tail surface 
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(Section 4.4 includes a discussion on the impingement Knudsen number and how 
it is used to determine the impingement regime). If the impingement mechanism Is 

I 

free molecular or transitional, the impinging molecules will be re-emitted from 
the surface in a random manner; some of them could impact the payload. If the 
impingement mechanism is continuum in nature, a boundary layer will be formed 
which will build up toward the top of the tail. Impingement regimes on the tail 
were calculated and the boundary layer height and mass flow in the boundary 
layer were calculated. 

4. 1.2. 4 Analysis Results 

The surface Knudsen number (which varies from 0.01 close to the tail base to 
0.0007 at the. tail tip) reveals the impingement mechanism of the aft +Z RCS 
plume on the tail is continuum in nature over the entire tall surface. Therefore, 
no molecules will be re-emitted from the surface. A boundary layer is formed 
(Figure 4. 1.2-1) which originates close to the tall base and builds as the plume 
flow moves up the tail surface. The maximum boundary layer thickness (which 
exists at the tail tip) is on the order of 17 inches; however, the mass flow in 
the boundary layer is on the order of 0.2 percent of the total mass flow of the 
RCS engine. 

4. 1.2. 5 Preliminary Conclusions and Recommendations 

Based on the above analysis, our conclusions and recommendations are: 

• Since the RCS plume impingement on the tail is continuum in nature 
and the mass flow in the boundary layer is such a small percentage of 
the total RCS mass flow, the effect of plume reflections on the payloads 
is negligible. Therefore, we recommend that reflection of RCS 
plumes not be included in the model. 
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FIGURE 4.1 .2-1. SKETCH OF CONTINUUM FLOW AROUND TAIL DUE TO 
AFT +Z RCS ENGINE PLUME 


4.2 Object Geometry Modelling Approaches 

The studies documented in this section deal with alternate formulations 
of the Object Geometry Module, OBGEOM. The following functions are 
allocated to OBGEOM (Reference 1): 

9 Decompose object into subshapes (e.g., cylinders, flat plates). 

• Decompose subshapes into "elemental flat plates". 

• Designate "contamination control points" for computation of 
contaminant flux. 

• Transform from subshape axes to thruster axes and from thruster 
axes to object axes. 

• Compute, if necessary, object shadowing (the shadowing of one 
part of the payload by another). 

• Spatially integrate the elemental forces and moments. 

• Compensate for any offset of the object moment center (e.g., 
c.g, offset). 

The following subsections describe studies which were conducted to answer 
the following questions: 

• How much can the subshape geometry representation be 
simplified (e.g., 12-sided polygon -^cylinder s^flat plate)? 

• How fine an elemental subdivision ("mesh size") is required for 
the model? Should the mesh size be fixed or variable? 

■ Is the capability to simulate object shadowing required in the 
model ? 

Basic object geometry data used in these studies are defined in 
Section 3.2.2. 
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4.2.1 Representation of Object Geometry 

4. 2. 1.1 The Question 

Is a drastic simplification of payload geometry representation acceptable 
in the new plume impingement model? Currently, there are three payload 
configurations under consideration: LDEF, MMS/SMM, and Skylab. Among 
these payloads, we only have detailed information on LDEF and Skylab 
geometries. Therefore, the MMS/SMM payload will not be -discussed in the 
present study. c 

4. 2.1. 2 Tradeoff Factors/Consideratlons, 

Implementation : In general, for the purposes of analysis, any complicated 

payload geometry may be represented and approximated' with varying levels of 
fidelity. The first step in the approximation would -be to replace the actual 
payload geometry by a collection of simple subshapes, e.g. , cylinder, ■ 
circular p'late, rectangular piate, etc. The next step would be to further 
approximate these simple subshapes by either performing a direct subdivision 
of each subshape into smaller elements, or making further simplifications -to 
the subshapes and then subdividing “he resulting shapes into elements. . 

Since detail information on the MMS/SMM configuration is unavailable at 
present, only LDEF and Skylab configurations -are presented in Figure 4. 2. 1-1. 

For the present study the actual LDEF configuration is either represented as 
a cylinder or a "T-plate" approximation. The Skylab configuration is 
represented either by a collection of circular -cylinders and rectangular plates., 
or approximated by boxes and rectangular plates as shown (see References 7 ‘ 
and 9). 

"Cost 11 Factors : The direct benefits of simplification of object representation 

are: fairly minor savings of .storage requirements and computer time, and a 

simplification of programming. There is also an indirect saving of computer 
time, due to the use of coarser mesh sizes (see Section 4.2.2) ‘appropriate 
with simplified geometries. 
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Actual LDEF- 



Solar panels 
seen edge-on 


-Approximation- by Cylinder 


"T-Plate" Approximation 


(a) LDEF CONFIGURATION __ 




(b) SKYLAB CONFIGURATION 


FIGURE 4. 2. 1-1 LDEF AND SKYLAB CONFIGURATIONS AND THEIR SIMPLIFYING REPRESENTATIONS 



Performance Factors : In general, a closer representation of an original 

payload geometry will give better accuracy in impingement calculations. On' 
the other hand, a drastic simplification of the payload geometry may induce 
substantial errors in the impingement calculations. 

4. 2. 1.3 Analysis Approach 

In order to study the effect of subshape representation, the PLIMP program • 
(Section 3.2.4) was used to calculate the disturbances for the LDEF geometry, 
approximated by either a circular cylinder or a T-plate configuration, at 
several locations and orientations. The impingement pressure used in PLIMP 
was the modified newtonian formulation. A finer mesh size was used for the 
cylinder representation than for the T-plate representation. 

In a related study (Reference 10), the Skylab configuration approximated by ■ 
a collection of cylinders and plates was analyzed by the PLIMP program; 
resulting forces and moments were compared for the same geometry approximated 
by boxes and plates in the QLDP program (Reference 7 ) , as reported in 
Reference 9 . 

4. 2. 1.4 Analysis Results 


Table 4. 2. 1-1 summarizes performance data for impingement calculations by 
PLIMP for the LDEF geometry. Cost data (not shown) indicated that running 
the T-Plate approximation on PLIMP takes roughly the same execution time as 
a cylinder representation with comparable mesh fineness. In a simpler pro- 
gram (where geometry calculations would be relatively more significant), the 
T-plate approximation might take a quarter to a half of the execution time 
for the cylinder. However, we note from Table 4. 2. 1-1 that errors are gen- 
erally large in both force and moment data for the T-plate approximation of 
the LDEF geometry. 

For convenience, we also repeat in Table 4, 2, 1-2 relevant results for Skylab 
which were reported in Reference 10. Again, large differences in force and 
moment data are noted between the two representations for the Skylab configu- 
ration. 
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TABLE 4.2. 1-1 PERFORMANCE VARIABLES FOR LDEF -REPRESENTATION 
BY CIRCULAR CYLINDER AND T-PLATES 


TEST CASE 
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3 
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-1.222 
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4 
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0.81 

20,31 

• 5 
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0.04 
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6 
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0. 
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NOTES: 

3 With respect to the cylinder results 





TABLE 4.2. 1-2 PERFORMANCE DATA COMPARISON OF TOTAL SKYLAB DISTURBANCE DUE' TO THREE 
UP-FIRING JETS: QLDP VS. PLIMP (from Reference 10) 


PROGRAM 

RAW PERFORMANCE DATA ' 

SENSITIVITY 3 
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(lb) 
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Force 

Moment 

AMag 

A- Dir 

■A -Mag 

Dir 

QLDP 

-25.69 

28.89 

- 

38.66 

II 

- 

wE 

603.27 

+64.9,^ 

3.65° 

hi 13. 8 % 

13.97° 

PLIMP 

-16.49 

16.64 

-0.68 

23.44 

o 

-68.11 

2S5E 

282.12 

- 

- 

- 



NOTES: 

a with respect to PLIMP results 












4. 2. 1.5 Preliminary Conclusions and Recommendations 


Based on the above analysis and results, we conclude’ that, although a 
certain degree of simplification' of payload configuration is required, a 
drastic simplification of geometry can induce significant errors in the 
impingement calculations,- Therefore, a drastic simplification of subshape 
representation is not recommended in the plume impingement model. 
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4.2.2 Element Subdivision Mesh Size 


4. 2. 2.1 The Question 

How fine an element subdivision "mesh size' 1 is required for ac- 
ceptable accuracy? Should the mesh size be fixed or variable? 

4. 2. 2. 2 Tradeoff Factors/Considerations 

Implementation : In the present study three different subshapes 

(cylinder, circular plate and rectangular plate) are considered. 
Each subshape is divided into a number of elements for impinge- 
ment force and moment calculations. Figure 4. 2. 2-1 shows how 
various subshapes are divided, using a single "fineness parameter", 
n. (The PLIMP program (Section 3.2.4) actually allows for two 
fineness parameters, n and m, for each subshape.) The fineness 
parameter(s) would be provided by user input. On-line varia- 
tion of mesh-size would be implemented by simple logic based on 
distance and/or angle "breakpoints". For each individual element, 
the program "quantizes" flowfield properties and impingement 
angle at the centroid (properties and angle are assumed uniform 
over the "elemental flat plate"). 

“Cost" Factors: Variation of n causes minor impact on computer 

storage and complexity of the programming. Execution time var- 
ies roughly as the square of n. 

Performance Factors : Using a finer 'mesh size should increase 

accuracy, by allowing more precise accounting for flowfield 
gradients and variation of impingement angle over the object 
surface. Figure 4. 2. 2-2 shows, qualitatively, the form of the 
expected cost/ performance effects of mesh-size variation for 
fixed object geometry and orientation, neglecting roundoff: 

• .At a fixed distance, increasing fineness will reduce 

error. 

• Increasing distance will reduce the fineness required 
to attain a given error level. 
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Cylinder 




Rectangular Plate 


FIGURE 4. 2. 2-1 SUBDIVISION OF VARIOUS SHAPES; 

CROSS INDICATES ELEMENT CENTROID 
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FIGURE 4. 2. 2-2 EXPECTED FORM OF COST/ PERFORMANCE 
VARIATION (FIXED GEOMETRY AND 
ORIENTATION; ROUNDOFF NEGLECTED) 
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4. 2. 2. 3 Analysis Approach 

To evaluate effects of subdivision mesh size, a number of runs 
of PL IMP were made, for both MMS/SMM and LDEF payloads (Sec- 
tion 3.2.2)at various on-centerline distances and orientations. 
Object shadowing computations (Section 4;2.3) were included 
for all data runs. Impingement pressure formulations used 
(Section 4.4.1) were modified newtonian-and oblique shock. 

In addition, to investigate the property gradient effect on the 
mesh size, a special case was also included with the LDEF being 
displaced 30 degrees (57.74 ft) off the nozzle center line at 
d = 100 ft. This position was chosen because the rate of change 
of dynamic pressure with angle is maximal around 30 degrees. 

Finally, three runs were made with a cylindrical payload 60 feet 
in length and 15 feet in diameter, which is the maximum payload 
size specified for the Orbiter. This cylindrical payload was placed 
in the RCS plume 30 degrees from the centerline. 

Three different mesh sizes were chosen for the cylinder surface 
and -its circular end plate, such that all the element dimensions 
are comparable for a particular run. The modified newtonian 
impingment formulation was used' for these three runs. 


4. 2. 2. 4 Analysis Results 

Table 4. 2. 2-1 provides case-by-case cost and performance data for 
the LDEF and MMS/SMM runs. Errors were surprisingly small in many 
cases. Execution time displays the expected trend. Data for the 
60 X 15 payload, shown in Table 4. 2. 2-2, also show little variation 
in accuracy as mesh' fineness is reduced. 
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TABLE 4. 2.2-1 COST i PERFORMANCE VARIABLES FOR 
ELEMENT SUBDIVISION NESH SIZE 


TEST CASE 

J2SBS1 

COST DATA 

RAH PERFORMANCE DATA 

SENSITIVITY 11 

Jtuober* 

d 

(ft) 

A 

(deg) 

0 

(deg) 

nxn 

Execution 

Time 

(Sec) 

Time 

R4ttO d 

F x 

(15) 

1 

■ K 4 

<lfr-ft) 

Force 

Errors 

Moment 

Error 

(i» 

A Hag 
(±S> 

A Dir 
(+.deg) 

HI 

m 


45 b 

60 x 60 

364.8 

- 

74.06 

46.62 

(ETTTM 

n 

_ 






30 x 30 

89.8 

0.25 

74 23 

46.50 

KL-lR.fl 

■Esm 

0.13 

0.04 





IS x 15 

25.9 

0.07 

MM 

mrxrm 

| -493.9 

0.26 

0 23 

0.18 





8x8 

. 7.2 

0.02 

74.49 

mrm 


0.18 

0.36 

0.12 





4x4 

2.8 

0 01 

70.95 

49.36 

-515 4 

1.23 

2.64 

4.16 

K 2 

SI 

m 

20 c 

30 X 30 

303.2 

r 

65.44 

21.73 

-117.8 

_ 

_ 

- 





ItW 

88 1 

0 29 

65.27 

Bfflj 

■noon 

WKtWSMM 

0.03 

1.44 





8x8 

22,0 

0,07 

65,34 

21,63 

nm 

0.17 __ 

0,05 

HKI^I 





4x4 

7.4 

0,02 

65.93 



mwrmt 

0.14 


M3 

m 

0 

70 c 


287.5 

- 

21.34 

2.74 

-71.1 

- 

* 

mam 





15 x 15 

63.6 

0.22 

20.99 



1.53 

0.49 ‘ 

0.42 | 





8x8 

22.2 

0.08 

19.65 

3.28 

-70.9 

■am 

2.16 






4x4 

6 8 

0.02 

19.08 

3.47 

-71.9 

10.13 

2.99 

1.12 

K4 

m 

0 

20 c 

30 x 30 

277.0 


16.71 

5.52 

-17.1 

- 

- 

- 





15 x 15 

64.3 

0.23 

16.68 

5.52 

-17.2 

0.17 

0.03 

KE9H 





8x8 

23.4 

0.08 

16.59 

6.52 

-17.6 

0.68 

0.12 






4x4 

6 6 

0.02 

16.56 

5 55 ' 

-17.9 

0.74 

0.25 

4.67 

LI 

im 

0 

0 b 

30 x 30 

35.2 

- 

BEE9H 

0.0 

'0.0 

- 

- 

- 






12.8 

0.36 


0.0 

0.0 

0.0 

0.0 

- 





8x8 

4.8 

0.14 

- 21.35 

0.0 

0,0 

0.09 

0,0 

- 





4x4 

1.9 

0.05 

21.41 

0.0 

0.0 

0.18 

0.0 

- 

' 12. 

iTTT 

0 

20° 

mi 

. 56.4 

_ 

18.81 

2.84 

6.08 


_ 

. 





15 x 15 

16.2 

0.29 

18.81 

2.84 

6.03 

0.0 

0.0 

0.82 





8x8 

4.7 

0.06 

18.85 

2.87 

5.85 

0.26 

0.07 „ 

3.78 





4x4 

2.1 

0,04 

13.70 

2.90 

£J3_ 

... 0.53 

0.23 

1 7.40 

13 

El 

0 

70 b 

im 

40.3 


27.71 

-8.80 

29.44 



_ 





15 x 15 

10 3 

wmm 

ZLJ& 

-8.79. 

PQ.Tt 

0-07 

0.01 

0.37. 





8x8 

3.6 

0.09 


-8.76 

23.54 

0.38 

0.01 

0.34 





4x4 

1.8 

0.05 

27.60 

-8.76 

27,24 

0.33 

0.01 

7.47 

L4 

W. 

0 

SO 6 

3Q x 30 

33.4 

- 

32.47 

0.0 

0.0 

- 

- 

- 





EXES 

8.8 

0.26 

32.55 

0.0 

0.0 

0.20 

0.0 

- 





8x8 

3.4 

0.10 

»>*!:■ 

0.0 

0.0 

0,30 

0.0 

- 





4x4 

1.7 

0.05 

32.20 

0.0 

0.0 

0.80 


- 

L5 

OI 


50° 

30 x 30 

33.1 

• - 


2.16 

SEES 

MEHii 

HEHH 

ml 





IS x IS 

9.23 

0.28 

9.114 

2.12 

27.53 

0.21 

0.26 

0.4 





8x8 

3.33 

0.101 

9.08 

2 12 

27 13 

0.56 

0.18 

1.9 





4x4 

1.80 

0.054 

7.57 

2.12 

20.77 

19.3 

2.33 

24.9 


NOTES: 

?H denotes WS/SMH, L denotes LDEF 
IHodffied newton fan Ipipin-geoent pressure 
^Oblique shock Impingement pressure 
■1th respect to the finest nesh size 


CEOHETftY* 
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TABLE 4. 2. 2-2 COST & PERFORMANCE VARIABLES FOR A LARGE PAYLOAD 
DUE TO ELEMENT SUBDIVISION MESH SIZE' " 


TEST CASE 

MESH FINENESS 

COST DATA 

RAW PERFORMANCE DATA 

SENSITIVITY 5 

Number 

d ' 
(ft) 

, ♦ 
(deg) 

0 

(deg) 

Cylinder 

nxnr 

Disc 

nxm 

Executi on 
Time 
(sec) 

Time. 

Ratio 0 

F x 

Ob) 

*y 

(lb) 

H z 

(Ib-ft) 

Force Errors 

Moment 

Error 

(«) 

A Mag 
(±5!) 

A Dir 
(±deg) 

Tl 

100 

30 

90 

30x60 

30x7 

79.40 

- 

21,531 

1.762 

239.20 

- 

- 

- 

T2 

100 

30 

90 

15x30 

15x4 

21.04 

0.265 

21.535 

1.762 

239.09 

0.02 

0.0009 

0.046 

T3 

100 

30 

90 

5x10 

5x2 

3.57 

0.045 

21.435 

1.759 

Bsnm 

0.44 

0.013 

1.96 


NOTES ; 

a n denotes angular subdivision; 
m denotes linear subdivision 

b with respect to the finest mesh size 


GEOMETRY; y 



x 



4. 2. 2. 5 Preliminary Conclusions and Recommendations 
Based upon the above analysis and results, it appears clear 
that there is no need to use a fine subdivision mesh. Some- 
thing of the order of 10 x 10 appears adequate in the "near 
field", and probably 5 x 5 in the "far field". It appears 
that the flow field property gradients do not contribute a 
significant effect in choosing mesh size. 

At a slight increase in program complexity, the mesh size 
can be made variable with payload distance and/or orientation. 
This feature is recommended for use with offline programs, 
where averag e computing time is important; it is of minor value 
in realtime simulations, where maximum computing time is most 
important. 



4,2.3 Object Shadowing 


4.2.3. 1 The Question 

Do we require an object-shadowing capability in the plume imping- 
ment model? 

Figure 4. 2.3-1 conceptually shows the types of object shadowing to 
be considered: front/back and subshape/subshape. At least half 
of every subshape is shadowed, since there is a "back side 11 . Fur- 
ther shadowing can occur with payloads which have multiple subshapes; 
one subshape may block the "line-of-sight" from another subshape 
to the thruster. This "optical" shadowing is the highest fidelity 
Implementation considered; no attempt is made to model the actual 
flow around an object. If shadowing tests are not implemented in 
the model, the flow is assumed to impinge on every element, regard- 
less of position and orientation. 

4. 2. 3. 2 Tradeoff Factors/Considerations 

To qualitatively assess "cost" and performance factors, first con- 
sider how object shadowing would be implemented in the model. 

Implementati on : Each subshape is divided into elements as depicted 
in Figure 4. 2. 3-1 (a); see Section 4.2.2. The outward normal of 
each element, rn is dotted onto the vector alonq the thruster line 
of sight to that element, If i • n_ <o, the element is visible; 
if H >o, the element is on the back side and is invisible to 
the plume. 

The subshape/subshape shadowing determination is more complex, 
as can be seen in Figure 4. 2. 3-2. The general methodology for 
determining shadowing of a particular element is as follows: 

1. Transform coordinates to write the equations of 
the subshape and line-of-sight in same axes. 

2. Determine whether an intersection of the subshape 
and line-of-sight exists within the limits of the 
subshape. 




(a) Front/back Shadowing 



(b) Subshape/Subshape "Optical" Shadowing 


Figure 4.2. 3-1 Types of Object Shadowing Considered 
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3. If there is an intersection, determine whether it is 
"upstream" of the current element. 

These calculations must be done for every element of every subshape. 
{Note that a solid cylinder is represented as three subshapes: a 
hollow cylindrical "tube" and two circular end plates.) Short 
cuts may, however, be possible. 

"Cost" Factors : The cost of the„ back side test is minor, since 

the element position must always be expressed in thruster coordi- 
nates to get local flow properties. The "cost" of subshape/sub- 
shape shadowing depends partially upon logical complexity (different 
form of equations for each type of subshape - e.g., cylinder, 
circular and rectangular plates). There are also core requirements 
for the additional equations and logic. The execution -time cost 
may be positive or negative, since computations of intersection 
equations (some of which may be quadratic) and coordinate trans- 
formations must be added, but flow parameters, forces and moments 
will not have to be determined for shadowed elements. 

Performance Factors : Neglecting back side shadowing would obviously 
cause unacceptable errors in the force and moment calculations. 
Neglecting subshape/subshape shadowing will also obviously intro- 
duce error into the force and moment computations— in general, in 
both magnitude and direction. One can easily conceive of "path- 
ological" cases for which neglecting subshape/subshape shadowing 
reverses the sign bf the computed moment; see Figure 4. 2. 3-3. 

4. 2. 3. 3 Analysis Approach 

The following steps were taken to evaluate tradeoff factors for 
object shadowing: PUMP runs (Section 3.2.4) were made at various 

on-centerline distances and orientations with respect to a single 
thruster, using MMS-SMM geometry data as defined in Section 3.2.2; 
see Figure 4. 2. 3-4 for a definition of coordinates. Each run was 
made with and without the subshape/subshape shadowing option enabled 
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/ 


■shadow**^ 


Figure 4, 2. 3-3 "Pathological 11 Case Showing Sign Reversal of Moment 



Figure 4. 2.3-4 Coordinates Illustration 
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(the front/back computation is not optional). Forces, moments, 
and execution time with and without shadowing were compared for 
each test case. The inviscid-flow flowfield used for all runs was 
provided by MOC (Section 3.2.3). All runs were made with a con- 
stant mesh size of 30x30 (see Section 4.2.2); the impingement form- 
ulation was sometimes modified newtonian and sometimes tangent 
wedge/oblique shock (see Section 4.4.1). 

4. 2. 3. 4 Analysis Results 

No comparative results of the effect of back side shadowing are 
provided, since this is not an optional feature in PUMP. 

Table 4. 2.3-1 provides case-by-case performance and ''cost' 1 data 
for the runs made to analyze subshape/subshape shadowing. Cost 
data is in terms of Univac 1108 Central Arithmetic Unit (CAU) 
execution time. 

The "fixed cost" (additional storage required) to implement sub- 
shape/subshape shadowing has been roughly estimated by identifying 
those PUMP routines which are used only for this computation, add- 
ing their core storage and converting- to decimal. The estimated 
total storage required for subshape/subshape shadowing In PUMP is 
approximately 800 1108 words. (Our implementation may not be the 
same as that in PUMP.) 

4. 2. 3. 5 Preliminary Conclusions and Recommendations 

Based upon the above analysis we have come to the following con- 
clusions : 

• Significant errors in force and moment determination 
may be expected if shadowing capability is not included 
in the model. 

• Shadowing routines will require additional core, 
but will -actually save computer time at some payload 
orientations. 

Therefore, we recommend that object shadowing, both back side and 
subshape/subshape, be implemented in the plume impingement model. 
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TABLE 4. 2. 3-1 COST AND PERFORMANCE VARIABLES FOR SUBSHAPE/SUBSHAPE SHADOWING 


j= 

<*D 


| TEST CASE 

"COST" 

DATA 


RAW PERFORMANCE DATA 

SENSITIVITY 

c 

No 

d 

(ft) 

9 

(deg) 

Execution Time 
(sec) 

Time 

Ratio 

/With) 

Without Shadowino 

WBRIEWSMBI 

Force Errors 

Moment 

Error 

(+*) 

F x 

Ob) 

F y 

Ob) 

M z 

(lb-ft) 

F x 

(1b) 

F y 

Ob) 

M Z 

(lb-ft) 

Without 

Shadowing 

With 

Shadowing 

AMag 

(±S) • 

ADir 
(1 deg) 

v Wo / 

1 

35 

wm 

101.1 

89.8 

0.80 

93.17 

50.20 

-414.8 

74.23 

46.50 

iEBH 

20.8 



2 

50 

isa 

63.3 

58.7 

0.93 

44.57 

23.01 

-152.6 

34,72 

21.77 

-184.5 

22.4 

4.78 

17.29 

3 

50 

mam 

400.4 

287.5 

0.72 

30.07 

1.93 

-69.0 

21.34 

2.74 

-71.1 

40.0 

3.64 

2.95 

4 

100 

20 a 

78.9 

79.5 

1.01 

17.55 

5.31 

-14.9 

16.77 

5.36 

-16.8 

4.1 

0.89 

11.31 

5 

100 

20 b 

323.8 

277.0 

0.86 

17.46 

5.54 , 

-14.7 

16,71 

5.52 

-17.1 

4.1 

0.68 

14.04 

6' 

100 

45 a 

69.2 

65.3 

0.94 

10.78 

5.23 

, -20.3 

8.14 

5.27 

-25.2 

23.6 

7.04 

19, ,44 

7 

100 

70 a 

71.5 

69.6 

0.97 

7.41 

0.77 

-11. 7- 

5.50 

0.54 

-16.3 

34.8 

0.33 

28.22 

’ 8 

100 

70 b 

427.6 

314.6 

0.74 

7.19 

0.21 

-11.1 

5.35 

0.14 

-14.2 

34.4 

2.92 

21.83 


Notes : 


a Modified Newtonian Impingement Pressure 
b Oblique Shock Impingement Pressure 
c Shadowing is considered the correct result 



Geometry: 









































4.3 Plume Flow-field Modeling Approaches 

This section deals with alternative formulations of the RCS engine plume flow- 
field module, FLOFLD. The FLOFLD module provides the plume thermodynamic and 
flow properties necessary to calculate the plume impingement pressure effects 
on the payload. The data provided by FLOFLD are derived from a standard thermo- 
chemical/method-of-characteristics generated plume flowfield which is discussed 
below. 

Since the Orbiter RCS engines are operational at orbital altitudes (above 100 
nautical miles), the plume will be highly expanded. As the plume expands, the 
density decreases, and therefore the intermolecular collisions decrease and the 
molecular mean free path increases. Eventually, the intermolecular collisions 
become so few that the gas no longer obeys the laws of continuum fluid mechanics 
during the plume expansion process. At the molecular level, a gas has three 
energy modes: translation, rotation, and vibration. In the continuum regime, 
equilibrium is maintained among these modes; however, this may not be true in 
the low density region of the plume. When the 'gas undergoes a sudden change 
in translational temperature, several collisions are required to bring the 
rotational temperature into equilibrium, and several thousand collisions are 
required to bring the vibrational temperature into equilibrium. If the. number 
of collisions is' not great enough to keep these energy modes in equilibrium; 
the law of • equipartition of energy no longer holds. Therefore, each energy 
mode "freezes" at a different location. The temperature associated with that 
particular energy mode is then frozen, and the specific heat ratio and gas 
constant are adjusted to account for the loss of the degrees of freedom asso- 
ciated with that particular energy mode. Once all modes of energy have frozen, 
the flow is considered free molecular. 

When the continuum flow is source-like, i.e., high Mach numbers (low Mach 
angles) the streamlines are nearly straight. If transition to noncontinuum 
flow begins in the source- like flow regime, the freezing of each energy mode 
at its particular location has a less pronounced effect on the plume flowfield. 
Therefore, in the plume definition presented here, only the transition to free 
molecular flow was defined, i.e., all energy modes have frozen. The flow prior 
to transition was assumed continuum, and the thermodynamic and flow properties 
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were defined using continuum gas dynamics. After transition to free molecular 
flow, free molecular flow techniques were used to characterize the plume flow- 
field. 


The degree of rarefaction in the plume flow depends on the ratio of the 
molecular mean free path in the plume, and a characteristic dimension 
used to characterize the transition mechanism. The characteristic dimension 
used here was chosen as the nozzle exit radius. The dimensionless ratio is 
the Knudsen number in the plume flowfield, K n p. The plume Knudsen number can 
be expressed in terms of familiar parameters of fluid mechanics. From the 
kinetic theory of gases, the coefficient of viscosity is related to the mean 
free path (References 11 and 12 ): 

* 0/2)\pJ n (4.3-1) 


where 

= plume coefficient of viscosity 
A^ = plume mean free path 
p M = plume density 

V m = plume average molecular velocity 

The average plume molecular velocity can be related to the plume speed of 
sound as follows. The speed of sound is given by: 

»» “ / Ir. PJ/P. ‘ (4-3-2) 

where P is the plume static pressure and is the plume specific heat 
ratio. The plume root-mean-square velocity of the molecules is given by 
(Reference 12) : 


V rms = * /3P 7P: 


(4.3-3) 


51 



The relation between root-mean-square velocity and average velocity 
(References 12 and 13) is: 


V m = yB/3 it ' V„ mc (4.3-4) 

m rms ' ' 

Confining (4.3-2), (4.3-3), and 4.3-4): 

V ra = /B/Vir ' a (4.3-5) 

Substituting (4.3-5) into (4.3-1), dividing by Lp (nozzle exit radius) and 
multiplying numerator and denominator by the plume flow velocity U^: 


X fpp M 

n P = LpJ-ir ^ 

where 

R = P °° U °° L P and ^ 

e u a„ 

^oo co 


(4.3-6) 


This is the plume Knudsen number in terms of familiar parameters of fluid 
mechanics. 

The calculation of the coefficient of viscosity in a reacting multispecies 
flow system, such, as exists in an exhaust plume, is quite complex. The 
coefficient of viscosity for species "i", (Reference 14) is: 


|i. « 2.67 X 10" 5 


*1 a < 2 > 2 >* 


where 

M,. = molecular weight of species "1" 
T ro = pi ume temperature 
a ^ = collision diameter of species "i" 

n (2,2)* _ j 3 -j na) »y collision integral 


(4.3-7) 
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For a mixture of gases with “j" components: 


I 

i=l 


1 + 


a ik 


k = 1 
k f i 


X i 


where 

= mole fraction of species "1" 
and G.^ is given by: 


G ,- k = 


[l + ^i/|i k ' (H k /Hj) 1/4 ]' 


2.83 /r+MT/f.^ 


(4.3-8) 


(4.3-9) 


The coefficient of viscosity for a reacting multispecies flow system can be 
calculated using the TRANS 72 computer program (Reference 15). This program 
was used to obtain the plume coefficient of viscosity used in Equation 4.3-6 
to define the plume Knudsen number. 

The criteria used here to define transition to free molecular flow in the 
plume is that suggested in Reference 11: 

^ > 3 (4.3-10) 

For the RCS plume, this indicates free molecular flow- is reached when the 
mean free path is approximately 4.4 times the nozzle exit radius. This is 
the value used to determine transition to free molecular flow in the RCS 
plume. 

Another important aspect which must be considered when defining the thermo- 
dynamic and flow properties in the plume is how the nozzle boundary layer 
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affects the plume flowfield. The reduced velocity in the boundary layer 
perturbs the plums expansion at the nozzle lip. This perturbation is quite 
acute at high altitudes, which is the case here. This effect on the plume 
expansion was accounted for using a turbulent boundary layer analysis 
starting at the nozzle throat, Reference 14. 

Real gas effects in the exhaust plume were accounted for assuming equilibrium 
chemical reactions throughout the flowfield utilizing a one-dimensional 
equilibrium chemical kinetics computer program (ODE) discussed in Section 
3.2.3. The thermodynamic properties and species concentrations for the RCS 
propellant system (MMH/^O^) are then used with a method-of-characteri sties 
(me) computer program (discussed in Section 3.2.2) to define the plume. 

The MOC is a mathematical technique which defines the plume flowfield along 
left-running characteristics, and stores the data in a cartesian coordinate 
system located at th'e nozzle exit. 

Figures 4.3-1 through 4.3-6 present Mach number, temperature, density, and 
dynamic pressure contours in the standard RCS engine plume both in the near 
field (up to 200 feet from the nozzle exit) and in the far field (up to 
1000 feet from the exit). The operating and geometric characteristics of 
the RCS engine are given in Section 3.2.1. 

The transition from continuum to free molecular flow and the reflected shock 
are also shown in the figures. After free molecular flow has been reached, 
the mass flow streamlines remain straight and the temperature, velocity, and 
liach numbers along these streamlines remain frozen at the values attained at 
the transition point. The labelling of each mass flow streamline represents 
the percent of the total mass flow in the plume which falls below that stream- 
line. Since the other thermodynamic and flow properties in the plume have 
frozen at the transition point and the streamlines are straight, the density 
simply changes as the inverse of the cross-sectional area (1/r law). 

This plume definition will be used to calculate the plume-induced pressure 
effects on the payloads. Since the plume flowfield data generated by the 
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method of characteristics are at random locations, an extensive search and 
interpolation scheme is required when flowfield information is needed at a 
particular location in the plume. It would therefore be very advantageous 
to be able to represent the plume using some type of closed-form technique 
which would make data retrieval fast and efficient. The possibility of 
using source flow equations or a mathematical curve fit technique to repre- 
sent the standard data is discussed in Section 4.3.1. 

'Many of the RCS engines whose plumes may affect the payloads have scarfed 
nozzles. The effect on the plume due to nozzle scarfing is discussed in 
Section 4.3.2. 

The engines may operate over a large range of altitudes (100 to 350 nautical 
miles). Effects of the ambient atmosphere in this altitude range are dis- 
cussed in Section 4.3.3. 

For certain maneuvers, engines may operate simultaneously. Plume charac- 
terization for simultaneous engine operation is discussed in Section 4.3.4. 
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.. Hach Number Contours In the RC5 Engine Plume (continuum and free molecular regimes) 
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FIGURE 4.J-2 • Temperature (*R) Contour* in the RCS Engine Plume (continuum and free molecular regimes) 
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FIGURE. 4.3-4 Density (slugs/ft ) Contours and Velocity (ft/sec) In the RC5 Engine Plume (free molecular regime) 
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FIGURE 4i3-5. Ojmjolc Pressure (lbf/ft ) Contours in the RCS fnqlno Plume (Continuum and free Macular regimes) 








4.3.1 Use of Source- Flow Approximations 


4. 3. 1.1 The Question 

What is the best method of representing the high-fidelity plume flowfield 
and obtaining data at desired points? 

The alternatives are: 

• Use some type of source flow equations to represent the MDAC-West standard 
plume data 

• Use the MDAC-West standard flowfield data described in Section 4. 3, either 
in raw tabular form, or reduced by conventional curve/surface-fit tech- 
niques. 

Various "source flow" and source flow-like models have been considered for 
representing plume flowfields. These models, described in some detail be- 
low, are as follows: 

• Simple source flow (entire flowfield) 

• Source flow with boundary-layer correction (entire flowfield) 

• Source-flow forms with "curve-fit" constants and corrections 
(entire flowfield) 

■ Patched source flow (flowfield subdivided into "regions") 

Simple Source How : The source flow model of a rocket plune (Reference 16) 
is derived from two basic assumptions: 

• At any. fixed angle from the thruster centerline, the gas acts as though 
emanating from a point source; thus, all streamlines are straight, and 
density along each streamline obeys an inverse-square law. 

• The thermodynamic parameter (ratio of specific heats) of the gas, y, is 
constant over the field. 

Choosing an empirical functional form for the angular dependence (Reference 
17), the dynamic pressure at any point (r,8) can then be calculated by: 

r iL cos 6 (c.0) 

Q r z (4.: 

where A, B, and C are constants determined by the thruster parameters and 
the choice of y. Similar equations can be developed for density, Mach 
number, and other flow properties. 
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Source Flow With Boundary Layer Correction : For a finite boundary layer 
thickness, the simple source flow analysis can be modified to include an 
exponential boundary layer correction term in the calculation of dynamic 
pressure at high angles from the plume centerline (Reference 17): 




CC6 B (C@) f 

<*os B (ce 0 )«* f bce-ejl 


(4. 3. 1-2) 


0>©o 


Here B and 0 O are additional constants which depend on the boundary layer 
thickness and nozzle exit radius. 


Source Flow With "Curve-Fit" Constants : This technique predicts dynamic 
pressures using the same form as Equation 4. 3. 1-2 *, however, the five con- 
stants are now considered mathematical "curve fitting" parameters which are 
calculated by curve fitting the standard data. 


MPAD Mathematical Curve Fit Model : MR AD has developed a similar method of 
calculating dynamic pressures based on curve fits to standard data. The 
dynamic pressure, q, is calculated as: 
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Patched Source Flow Equations : With this technique, we again use Equation 
4, 3, 1-2 with "curve-fit" parameters. However, the constants are determined 
independently in various "regions" of the plume, rather than for the entire 
plume at once. (Of course, fitting .the entire plume can be considered a 
special case of this procedure.) These regions are then "patched" together 
to cover the entire flowfield, generating a table of source-flow parameters 
as functions of thruster coordinates. 

Note that this technique may cause non-physical discontinuities in flowfield 
properties at the region boundaries. It may be desirable to use overlapping 
patches and/or some type of averaging scheme. 

Manipulation of Standard Data : MDAC-West standard data can be stored in 
on-line tables; search and interpolation techniques can then be used to 
find flowfield properties at given thruster coordinates. Another approach 
Is to fit the standard data with conventional mathematical curves or sur- 
faces— e.g. , polynomials or splines, 

4. 3. 1.2 Tradeoff Factors/Considerations 


Implementation : Implementation of the entire-field source flow and source 

flow-like techniques is simple and straightforward. Implementation of the 
patched source flow technique will require searching a small data block in 
the program. Using the technique involving manipulation of standard data 
will require complex search and interpolation schemes and logic in the program. 

In implementing a flowfield model for plume impingement calculations, one 
must include a streamline model in order to calculate impingement angles. 

For the high-fidelity streamlines shown in Section 4.3, the coirmon assump- 
tion of radial streamlines appears justified for all distances from the 
nozzle large compared to the nozzle exit diameter. 

"Cost" Factors : . In terms of complexity, storage, and execution time, the 
simple source flow, source flow with boundary layer correction, and source 
flow with curve fit constants are the easiest and fastest to use, since 



they require at most one comparison (8 vs. e 0 ) and one calculation. The 
MPAD model and the patched source flow are also simple, and have fast 
execution times. Direct use of the MDAC-W standard data will require 
very large storage areas, and considerable complexity in manipulation 
of data. Storage and execution time requirements for the use of con- 
ventional- curve fits to the standard data are difficult to estimate at 
this time. Execution time is expected to be proportional to the com- 
plexity, with one special caution; if the available main-storage area 
were too small to hold all the standard data, the need to use "overlay” 
techniques would increase execution time dramatically. 

Performance Factors : Since the expected end use of this model will include 
analysis of "R-bar" and "V-bar" approaches (Reference 18), the plume flow- 
field must be accurately modelled in the highly rarlfied regions of the 
plume at large angles from the centerline, even at the expense of storage 
and execution time. Accuracy requirements are of fundamental importance 
in recommending a plisne flowfield technique. Those flowfield models- 
which satisfy certain accuracy criteria will then be considered for min- 
imum storage and execution time. 

4.3.1. 3 Analysis Approach 

To calculate the theoretical constants in a- source-flow model, a value must 
.be chosen for y, all other parameters required are engine parameters (e.g., 
chamber pressure). For the simple source flow, y must be a constant; we use 
Y = 1.236, as suggested in Reference 19. Actually, y varies in the plume. 
(In regions of low chemical activity, y is basically a function of gas tem- 
perature.) Therefore, it is important to consider the sensitivity of the 
source-flow results to the value of y. 

The constants in Equation 4. 3. 1-1 are all functions of y. Computing these 
constants for various y's allows one to plot dynamic pressure contours 
using simple source flow analysis. Figure 4. 3. 1-1 shows two example dy- 
namic pressure contours for various y's. Qualitatively, we see that the 
value of y can significantly affect the pressure contours. 

To get more quantitative results, the radial dependence in Equation 4.3. 1-1 
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Legend: 

A y = 1.20 
B y = 1.22 
C y » 1.236 
D y = 1.26 


E y = 1.28 
F Y “ 1.30 
G Y “ 1.32 
H y c 1.34 



FIGURE 4. 3. 1-1 SENSITIVITY OF DYNAMIC PRESSURE CONTOURS TO VALUE OF THERMODYNAMIC PARAMETER ( Y ) 


2 

is eliminated by plotting qr as a function of y for various values of 
streamline angle, 0, as shown in Figure 4. 3. 1-2. Since the plots are rough- 
ly linear for a given 0, we see that the sensitivity of the dynamic pres- 
sure to y is constant for a given 8. On the centerline (8 = 0) a 1% change 
in y gives a 3% change in dynamic pressure; at 0 = 60°, a 1% change in y 
gives a 15% change in dynamic pressure. Therefore, the dynamic pressure 
contours generated using source flow analysis are very sensitive to y. 

The various methods for generating flowfields described in Section 4.3.1-.1 
were used to determine dynamic pressures. The results were compared to 
MDAC-W standard dynamic pressure data. 

Two sets of comparisons were made to the MDAC-W data. The first comparison 
is qualitative, and consists of overlaying dynamic pressure contours. The 
second comparison is a quantitative analysis of radial and angular slices 
of dynamic pressure. In investigating the patched source flow technique, 
a few representative sample regions were chosen for evaluation. 

4. 3. 1.4 Analysis Results 

Constants used in the five techniques are given in Table 4. 3. 1-1. Dynamic 
pressure contour overlays for the various techniques are shown in Figures 
4. 3. 1-3 through -7; in each figure, (a) is the near-field analysis, (b) is 
the far field. Our data indicate that a source flow model must include a 
boundary layer correction in order to predict high-angle pressures accurately. 

The comparisons to MDAC-W dynamic pressure data along radial and angular 
slices are shown in Figures 4. 3. 1-8 through 4.3.1-12. We believe the 
re-plotted MDAC-W data to be accurate within +2 feet. This random plotting 
error may be significant at high angles, where the pressure gradient is 
steep, (To generate actual model data, we would use the standard data 
directly, in digital form, thus eliminating such plotting errors.) 

The simple source flow model appears good at low angles, with typically 
10% error, but fails at high angles. The technique using source flow with 
a boundary layer correction gives better agreement (at least in the form of 
the contours) at large angles, but errors in this region. are still large. 
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FIGURE 4. 3. 1-2 SENSITIVITY OF SOURCE FLOW - FLOWFIELD TO yi RADIAL SLICES 


TABLE 4.3. 1-1 CONSTANTS USED 


' ^^CONSTANTS 
MODEL 

A 

B 

Simple Source Flow 

480.4 

8.47458 

Source Flow with 
Boundary Layer 
Correction 

480.4 

8.47458 

Source Flow with 
Curve-Fit Con- 
stants 

480.4 

8.47458 

MPAD Fit 

N/A - see Equation 4. 3. 1-3 

Patched Source 
Flow: 

Region a (1) 

500 


(2) 

480 


. (3) 

480 

8.47458 

(4) 

400 


‘ .(§L_ 

460 



a 


See Figure 4. 3. 1-7 































FIGURE 4.3. 1-3 SIMPLE SOURCE FLOW VS 
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(b) Far Field Dynamic Pressure Contours 


FIGURE 4. 3. 1-3 (Continued) 
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(a) Near Field Dynamic Pressure Contours 


Axial distance (ft) 


FIGURE 4. 3. 1-4 SOURCE FLOW WITH BOUNDARY LAYER CORRECTION VS. STANDARD DATA 
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(a) Near Field Dynamic Pressure Contours 


FIGURE 4. 3.1-5 SOURCE FLOW WITH CURVE FITTING CONSTANTS VS. STANDARD DATA 
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(b) Far Field Dynamic Pressure Contours- 


FIGURE 4. 3. 1-5 ‘ (Continued) 
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FIGURE 4.3. 1-6 MPAD FIT VS. STANDARD DATA 
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(b) Far Field Dynamic Pressure Contours 


FIGURE 4,3. 1-6 


(Continued) 




FIGURE 4.3. 1-7 PATCHED SOURCE FLOW VS. STANDARD DATA 
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(b) Far Field Dynamic Pressure Contours 


FIGURE 4.311-7 (continued) 
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(b) Angular Slices 


FIGURE 4.3. 1-8 (continued) 
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(b) Angular Slices 
• FIGURE 4. 3. 1-9 (continued) 
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dynamic pressure (psf) 



FIGURE 4.3.1-10 SOURCE FLOW WITH CURVE FITTING CONSTANTS VS. STANDARD DATA 
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(b) Angular Slices 


FIGURE 4.3.1-10 (continued) 
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dynamic pressure (psf) 



(b) Angular Slices 

FIGURE 4.3.1-11 (continued) 
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(b) Angular Slices - 
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FIGURE 4.3.1-12 (continued) 
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The source flow technique with curve-fit constants shows good agreement, 
with typical errors less than 10%. The MPAD model shows similar accuracy, 
with slightly larger core requirements. The MPAD model fails in the high 
angle, very close-in region, but this region is of little concern to Shuttle 
plume problems. 

The technique which uses patched source -flow- like equations typically has 
errors around 5%, except for high angle close-in distances. It is interest- 
ing to note that only the scaling factor. A, changed in the different re- 
gions of the plume under consideration. 

4.3. 1.5 Preliminary Conclusions and Recommendations 
Based upon the above analysis and results, techniques using theoretical 
values for source flow constants and boundary layer constants do not pre- 
dict dynamic pressures to within 10% of the standard data for all regions 
of the -plume, and dynamic pressures are seen to be highly sensitive 
to the chosen value of the constant y. We also conclude that a boundary 
layer correction term is necessary. 

The source flow technique with curve-fit constants and. the MPAD fit both 
predict standard data to about 10% accuracy. Using patched source flow 
equations, errors. are typically less than 10%. 

It appears that the patched source flow technique is the most promising as 
far as minimizing complexity and maximizing accuracy. Therefore, this 
method is recommended for the flowfield model. 

The assumption of radial streamlines is justified in the regions of the 
plume under consideration and is also recommended. 
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4.3.2 Thruster Scarfing (Short Side/Long Side) 


4. 3. 2.1 The Question 

Is it necessary to account for the scarfed nozzle effect on the plume 
definition? The engines which must be considered are the engines whose 
plumes (or partial plumes) may affect payload deployment and retrieval 
operation, i.e., the plumes which expand above the orbiter vehicle. These 
engines are*. 


Engine 

Plume to be Considered* 

Forward +Z 

Entire plume 

Forward -X 

Only the portion of the 
plume which expands from 
the short side of the 
scarfed nozzle 

Forward tY 

Only the portion of the plume 
which expands above the orbiter 

Aft tY 

Only the portion of the plume 
which expands above the orbiter 


*See Figure 4.3. 2-1 

The -X engine is the most severely scarfed of the engines of concern here; 
however, the region of the plume which will affect the payloads expands 
from the short side of the scarfed nozzle;' see Figure 4. 3. 2-1. The scarfing 
of the short side of the -X engine nozzle is quite minor in comparison to 
the standard unscarfed nozzle', see Figure 4. 3.2-2. The forward +Z engine 
nozzles are scarfed the most of the engines whose plumes or partial plumes 
are of concern here. However* the scarfing of these engines is also quite 
small; the extent of scarfing for these engines is also presented in 
Figure 4. 3. 2-2. An indication of the extent of nozzle scarfing for the 
forward +Z engines and the forward and aft tY engines can be realized by 
noting the thrust vector is altered no more than 1.52° for any of these 
engines. The thrust vector for the -X engines is altered about 5.6° due to 
nozzle scarfing; however, the majority of this is due to the long side of 
the scarfed nozzle which does not affect the plume of concern here. 
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FIGURE 4. 3.2-2 FORWARD +Z AND -X RCS NOZZLE SCARF CONFIGURATIONS 
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Before this question can be answered, an understanding of the effect on the. 
plume characterization due to nozzle scarfing is necessary. The scarfing 
of the nozzle alters the flow at the exit of the nozzle, and in turn the 
plume expansion at the nozzle lip. The flow will expand to a lower exit 
pressure at the long side of the scarf due to the increased expansion ratio, 
and in turn renders a smaller Prandtl -Meyer expansion at the nozzle lip. 

Also, the smaller lip angle on the long side renders a smaller total expan- 
sion. Therefore, the scarfing of these nozzles results in a reduced expan- 
sion at the nozzle lip. This reduced expansion at the nozzle lip affects 
only the outer regions of the plume; the internal regions of the plume are 
essentially unaffected by nozzle scarfing. 

4. 3. 2.2 Tradeoff Factor/ Considerations 

To assess cost and performance factors, consideration of how the plume flow 
field from a scarfed nozzle would be implemented in the model is necessary. 

Implementation : To implement the scarfed engine plumes in the model, a 

separate plume characterization for each scarfed engine must be generated 
and stored. This would require the storage of several plumes in addition 
to the unscarfed engine plume. Additional logic would be required to 
establish if the operational engine (or engines) are scarfed or unscarfed 
and then, since the scarfed engine plumes are not symmetrical, the orientation 
of the plume would need to be established. 

“Cost 11 Factors : The cost of implementing the effects of the scarfed nozzles 
is the cost associated with storing, selecting, and accessing the scarfed 
engine plume flow field definitions, and establishing plume orientation. 

Performance Factors : Since the effect on the plumes of interest here due to 
nozzle scarfing is minor, the increase in accuracy in defining the plume 
Impingement effects utilizing the scarfed engine plume as opposed to utili- 
zing the standard unscarfed engine plume is minimal. 

4. 3. 2.3 Analysis Approach 

The forward +Z engines have the most severely scarfed nozzles of the engines 
whose plumes are of concern here (the outside engines in this cluster being 
scarfed the most). Dynamic pressure contours were generated and compared 
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for the standard unscarfed .RCS engine plume and the outside forward +Z 
engine plume. Figure 4. 3.2-3. 

4. 3. 2. 4 Analysis Results 

The upper half of Figure 4. 3. 2-3 presents the dynamic pressure contours in 
the standard unscarfed RCS engine plume and the lower half of the plot 
presents the same information for the long side of the outside forward +Z 
scarfed engine. Comparison of these contours indicates .that the effect on 
the plume definition due to nozzle scarfing for these engines is minor, 
and in turn the effect on the plume impingement results due to nozzle scarfing 
is negligible. Since the nozzle scarfing for the lY engines (both forward 
and aft.) is less severe than the long side of the forward +Z engines, the 
effect on those plumes would be less pronounced than that presented in 
Figure 4. 3. 2-3. 

4. 3. 2. 5 Preliminary Conclusions and Recommendations 

Based on the above analysis, our conclusions and recommendations are: 

•The effects on the plume flow field definition due to nozzle scarfing 
for the plumes (or partial plumes) of concern for payload deployment 
and retrieval should be neglected. A single plume definition based on 
the standard unscarfed RCS engine should be used for all the RCS engine 
■plumes (scarfed and unscarfed) of interest here. 





4.3.3 Ambient Atmosphere Effects 
4.3. 3.1 The Question 

Does the plume impingement model require a capability to account for high 
altitude ambient atmosphere effects? There are two topics related to this 
question: first, the interaction between an exhaust plume and the ambient 

atmosphere; second, the interaction between the ambient atmosphere and a 
payload. 

Plume/atmosphere interaction : We are interested only in the orbital altitude 

range of 100 to 350 nautical miles (nm); 120 nm is considered "typical". 

At these high altitudes the ambient air is streaming past the Orbiter and 
payload at an orbital speed of approximately 25 kfps; however, its density 
is so rarefied that its molecular mean free path, A (which is inversely 
proportional to the air density, p) is very large. For example, at 120 nm, 
the air mean free 'path is greater than 1000 feet, as can be seen from 
Table 4.3. 3-1, taken from Reference 20. Thus, the molecules of the air 
experience very few collisions (i.e., the atmosphere is in free-molecule 
flow) at this highly rarefied condition. 

On the other hand, a rocket exhaust plume resulting from expanding into a 
vacuum environment may be considered as consisting of a continuum core, 
followed by a transition region and then by a large rarefied or free-molecule 
zone, as pictured in Figure 4. 3. 3-1. (Note that, in reality, the boundaries 
between these flowfield regimes are not sharp as shown; in fact, flowfield 
properties vary smoothly across these artificial boundaries.) 

When the rarefied ambient atmosphere and the vacuum plume described above 
encounter each other at high speed, we would expect some kind of interaction 
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TABLE 4.3. 3-1 ATMOSPHERE CHARACTERISTICS AND ORBITAL VELOCITY 
'FOR ALTITUDES OF INTEREST 


ALTITUDE 
(hj nm) 

AMBIENT ATMOSPHERE 
CHARACTERISTICS 5 

ORBITAL VELOCITY 
(ucoj ft/ sec) 

MEAN FREE PATH 
(X, ft) 

DENSITY „ 
(p» slug/ft) 

100 

475.1 

9.77 x 10' 13 

25043.0 

120 

1209.7 

3.68 x 10" 13 

24975.7 

150 

3875.0 

1.07 x IQ-" 13 

24875.5 

200 

18169.0 

2.04 x 10.-™ 

> 24711.0 

350 

561680.0 

5.24 x TO-™ 

24236. 5 


a From Reference 20 



FIGURE 4. 3. 3-1 SKETCH OF FLOWFIELD REGIMES FOR 
A ROCKET PLUME EXHAUSTING INTO 
A VACUUM 
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to occur. We wish to determine whether this interaction distorts the 
plume flow-field as compared to the "free-, space" (vacuum expansion) form 
assumed in generation of the "standard" data of Section 4.3. 

Atmosphere/ pay! oad interaction : When a payload is placed at close range 
to the vacuum plume, the dynamic pressure due to the plume is much greater 
than that of the ambient atmosphere. As the range is increased, the plume 
dynamic pressure falls off rapidly (roughly as p?; see Section 4.3.1). 
Eventually the plume dynamic pressure will become "negligible" as compared 
to that of the ambient air. 

Therefore, we could establish - a criterion for a "practical boundary" of 
the plume: say, q rcs < 0.1 q aepQ , where q rcs and c « aero denote plume 

and air dynamic pressures, respectively. We wish to determine whether 
this "practical boundary" falls in the range of interest for the plume 
impingement model (35-1000 feet) or not. 

4. 3.3. 2 Tradeoff Factors/Considerations' 

Implementation : If we establish that the plume distortion due to the ambient 

air is significant, modelling this effect will require offline computations 
of the plume distortion boundary, and then storing the boundary data in 
table/curve form, e.g., as functions of altitude, distance, and angle. 
Additional model inputs and logic would be required to establish the oper- 
ational altitude and select the appropriate plume definition. 

On the other hand., if we establish that the "practical boundary" from the 
dynamic pressure ratio falls in the range of interest, it can be implemented 
with a simple test, using parameters obtained from a fairly simple offline 
analysis. 
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"Cost" .Factors : Even assuming that we could find or develop means for 

the required offline computation of plume distortion, the associated fixed 
cost would be difficult to estimate. In addition, the online modelling 
cost could be substantial. 

On the other hand, implementing the "practical boundary" -costs little in 

t 

either offline or online modelling. It may actually save online computer 
time, since plume impingement computations are bypassed whenever the test 
fails. 

Performance Factors : Plume distortion, if it exists, is only on the 

"windward" side of the plume. For most operations of current interest 
{e.g., nominal braking, R-bar, and V-bar approaches), the payload will 
usually be on the "leeward" side of the plume. Therefore, there seems to. 
be little or no accuracy payoff from modelling plume distortion. 

4. 3. 3. 3 Analysis Approach 

Based upon the above considerations and the fact that obtaining a quantitative 
assessment of the high altitude effects on the plume flowfield would 
require a substantial effort involving an analysis of two highly rarefied 
interacting gases, we decided that only a qualitative assessment of the 
interaction was warranted. To obtain such an assessment, we were led to 
conduct this part of study by consultation with experts in the field and a 
survey of the literature. At least a half dozen aerodynami cists were 
contacted, and a dozen related papers were identified and reviewed. 

For the "practical boundary" we computed the ambient dynamic pressure, 
a t various altitudes, using the orbital speeds and atmospheric 

aciO 
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densities shown in Table 4. 3. 3-1. Then- we compared these dynamic pressures 
with the plume dynamic pressure generated by MDAC-West as reported in 
Section 4.3. 

4. 3. 3. 4 Analysis Results 

PI ume/atmosphere i nteracti on : None of the plume experts we consulted had 

actually studied the problem in detail, but most of them were confident 
that the interaction is negligible {from the viewpoint of the plume being 
affected by the air). The same conclusion has also been reached from 
the literature survey. 

The literature survey also uncovered an analysis of the effect of the plume 
upon the air at high altitudes. A rather complex technique available for 
computing this effect is reported in Reference 21; a physical- picture of 
the interaction will be qualitatively described below. 

With the vacuum plume sweeping at high orbital speed, the ambient air 
first encounters the outer portion of the plume. Since both the plume and 
air are highly rarefied, they pass through each other undergoing occasional 
collisions, without being disturbed much. As the ambient air penetrates 
further into the denser portion of the exhaust plume, the air stream is 
scattered by the exhaust gas, and some of the scattered air molecules will 
be deposited with and carried along by the gas. In the process the number 
density of the air is increased gradually as compared with the unscattered 
state, without affecting the plume. 

i 

Finally, as the air stream reaches the high density plume in the vicinity 
of the continuum core, the air molecules can no longer penetrate this region, 
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and are completely scattered by the denser gas. During the last stage of 
penetration the air number density peaks out and then drops off rapidly, 
while locally the plume density increases very slightly in the process. 

In an example calculated for a plume at roughly 350 nm altitude, it was shown 
that the ambient air can penetrate deep into the plume, and that its 
density may increase by up to a factor of ten at high angles from the plume 
centerline, due to the plume scattering effect. However, even a ten-fold 
increase in the air density does not affect the plume at all in this case. 

Atmosphere/payload interaction : Segments of the approximate dynamic 

pressure equivalence boundaries (dynamic pressure ratio q /q = 1, 

res aero 

using the MDAC-West dynamic pressure) for several altitudes are presented 
in Figure 4, 3. 3-2. Note that no plume effect on the air is accounted for 
in this plot. It seems clear that if a "practical boundary" -is defined 
at' a dynamic pressure ratio of 0.1, then this boundary will fall within 
the distance range of interest for altitudes of interest. 

4.3.3. 5 Preliminary Conclusions and Recommendations 

From the above analysis, we conclude that: (a) plume distortion by the 

S 

rarefied atmosphere is negligible, (b) the air dynamic pressure becomes 
comparable to the plume dynamic pressure for altitudes and distances of 
interest. 

Therefore, we do not recommend implementation of a plume distortion 
computation, but do recommend implementation of the "practical boundary" 
test in the plume impingement model. 
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4.3.4 Multiple Simultaneous Firings 
4.3.4. 1 The Question 

What is the best method of handling multiple simultaneous engine operation? 

This question must be answered for: (1) simultaneous operation of adjacent 

engines (same cluster) and (2) simultaneous operation of separated engines 
(different clusters). The word "cluster" is used here to indicate engines in 
the same pod which fire in the same direction. 

The alternatives are: 

• Considering each engine independently of the other and simply summing 
their effects. 

• Using a special technique to account for the interaction of the 
multiple plumes. 

Before this question can be answered, an understanding of the plume interaction 
structure is necessary. 

Simultaneous Operation of Adjacent Engines : Figure 4. 3. 4-1 depicts the flow 

structure created when two adjacent RCS engines are operated simultaneously. 

As the adjacent plumes intersect, a slip line (or plane) is formed midway 
between the engines. This plane acts as a solid boundary to the interacting 
plumes. Since the flow must turn to become parallel to the slip plane, an 
interaction shock wave is formed in each plume. This shock strength dissipates 
with increasing distance from the noz 2 le exit because the flow turns less, 
the greater the distance from the nozzle exit. In a relatively short distance 
(much less than the minimum distance from the engine exit to the affected- 
body), the interaction shocks approach the respective plume centerlines and 
dissipate completely. After the shocks have dissipated, the plume structure 
from the two adjacent plumes is more accurately represented by a plume from a 
single "equivalent" engine located midway between the single engines. The 
"equivalent" engine is a scaled up version of the single engine to render the 
same mass flow as the two single engines operating simultaneously. 
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SKETCH NOT TO SCALE 



FIGURE 4. 3.4-1 MULTIPLE ENGINE PLUME INTERACTION CHARACTERISTICS AND EQUIVALENT ENGINE PLUME 



Simultaneous Operation of Separated Engines : For the separated engines two 

potential plume interaction situations must be considered: (1) interaction 

between the up firing (+Z) engine plumes in the forward and aft pods and 
(2) interaction between the up firing (+Z) engine plumes in the left and 
right aft pods. The plumes from the (+Z) engines in the left arid right 
■aft pods will not interact because of vehicle geometry. The tail effectively 
blocks these plumes before they can interact. The plumes from the +Z engines 
in the forward and aft pods will interact. The interaction takes place at 
large enough distances from the engine exits that the plume flow fields are 
free molecular in nature; therefore, the existence of one plume will not 
affect the other. 

4. 3.4. 2 Tradeoff Factors/Considerations 

To assess tradeoff factors, the frequency of simultaneous engine operation 
must be considered. The frequency of simultaneous engine operation depends 
on the type and number of maneuvers requiring multiple engine operation, 
and is difficult to evaluate at this time. Therefore, frequency of multiple 
engine operation must be considered at a later date. To assess cost and 
performance factors, consideration of how simultaneous engine operation would 
be implemented "in the model is necessary. 

Implementation : Since the plumes resulting from simultaneous operation of 
two adjacent engines can be represented by a single “equivalent" engine, only 
a single plume flow field need be considered when calculating the multiple 
plume effects. The "equivalent" engine operating and geometric characteristics 
are identical to the single RCS engine except the "equivalent" engine is/2 
larger to obtain the same mass flow as the two single engines operating simul- 
taneously* This eliminates the necessity of generating and storing an addi- 
tional plume definition for simultaneous engine operation; a simple scaling 
of the single engine plume is all that is required. 

To implement the "equivalent" engine to simulate adjacent simultaneous oper- 
ation, additional logic would be required to determine whether the operational 
engines are adjacent or separated. If the engines are adjacent, then the 
"eouivalent" engine plumes would be used. If the engines are separated, then 
each plume would be considered independently and their effects summed. 

*The /2~renders the "equivalent" engine exit area equal to the sum of the two 
single engine exit areas. 


107 



A second alternative in implementing simultaneous adjacent engine operation 
is to consider each plume independently and sum the results. This technique 
would eliminate the requirement of additional logic to test whether the 
engines operational are adjacent. 

"Cost" Factors : If a single "equivalent" engine plume is used to represent 

the two adjacent engine plumes, the calculation time to establish the plume 
effects is of the order of half the time required to establish the plume 
effects from each adjacent engine. The additional logic necessary to imple- 
ment the "equivalent" engine simulation is minimal in comparison to the total 
plume effects computation time; therefore, substantial savings may be realized 
if simultaneous firing of adjacent engines occurs frequently. 

For the case of separated engines, no savings in computational time or storage 
costs can be realized since each plume and its effects must be considered 
separately. 

Performance Factors ; The plume flow field structure resulting from simul- 
taneous operation of adjacent engines is represented quite accurately in the 
region of interest here by the "equivalent" engine plume (see Figure 4. 3.4-1)'. 
However, since the spacing of the adjacent Space Shuttle RCS engines is only 
a few inches, the error induced in neglecting the multi pie plume interaction 
and simply surrming the results from each is minimal. 

4. 3. 4. 3 Analysis Approach 

Dynamic pressure contours and the free molecular freeze locations in a single 
engine plume and in an "equivalent" engine plume were generated and compared; 
see Figure 4. 3.4-2. The dynamic pressure contours are typical; other plume 
flow field variables (Mach number, pressure, density, etc.) can be compared 
in a like manner. 

4. 3. 4. 4 Analysis Results 

The sft scale factor relating the single engine plume and the "equivalent" 
engine plume can be verified from the dynamic pressure contours in Figure 4. 3. 4-2. 
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This scale factor is typical for all plume flow field variables (Mach number, 
pressure, density, etc.). 

4. 3. 4. 5 Preliminary Conclusions and Recommendations 

Based on the above analysis, our conclusions and -recommendations are: 

• If simultaneous adjacent engine operation occurs frequently, the 
calculation of the plume effects utilizing an "equivalent" engine 
is recommended. If simultaneous adjacent engine operation occurs 
very infrequently, then calculation of the plume effects from 
each engine independently of the other and summing the results is 
recommended. 

• During simultaneous operation of separated engines, the calculation 
of the plume effects from each' engine independently of the other 
and summing the results is recomnended. 




4.4 IMPINGEMENT PRESSURE MODELLING 


This section deals with alternative formulations of the plume impingement 
pressure module, IMPRES. The IMPRES module calculates the local normal and 
tangential plume impingement pressure coefficients, which are used to define 
the plume-induced forces and moments on the payloads. 


As described in Section 4.3, the exhaust plume at orbital altitudes will be 
highly expanded. In the far field of the plume, the flow will become highly 
rarefied; i.e., the fluid, does not behave as a continuous fluid (continuum gas 
dynamics) but rather exhibits characteristics of its molecular structure (rare- 
fied flow). 

When the low density portion of the plume impinges the surface of a payload, 
consideration of the ratio of the molecular mean free path in the plume, 
to a characteristic dimension of the impinged surface, L s , is necessary. This 
dimensionless ratio is the Knudsen number of the impinged surface, K^. 


The Knudsen number in terms of familiar parameters of fluid mechanics is derived 
in Section 4.3 and will not be repeated here. The Knudsen number used to define 
the degree of rarefi cation of the plume impingement mechanism on the surface is: 


where 





p U L e 
r «o « s 



and M = ^ 


(4.4-1 } 


1 1 1 



Rarefied flow techniques must be employed when the impingement Knudsen number 
becomes larger than about 0.01. The mechanism by which the impingement pressures 
are calculated is divided Into three impingement regimes depending on the degree 
of rarefaction (i.e. the impingement Knudsen number). The impingement regimes. 
In general, are defined as: 

0 £ K n$ < 0.01 continuum impingement 

0*01 £ K ns < 10.0 transitional impingement 

10.0 £ K ns < “ free molecular impingement 

Figure 4.4-1 presents these impingement regimes for a 10 foot payload surface 
(i.e. L = 10 ft.) in the RCS plume. Impingement regimes for any size payload 
surface in the RCS plume can be obtained from Figure 4.4-1 utilizing the defi- 
nition of the impingement Knudsen number and the impingement regimes defined 
above. 

Each Impingement mechanism is discussed separately In Sections 4.4.1, 4.4.3, 
and 4.4.2 respectively. Basically, in the continuum regime the impinging flow 
behaves as a continuous fluid; in the transitional regime the flow is only 
moderately rarefied, so that the molecules being reflected from the surface 
affect the oncoming molecules; and in the free molecular regime the flow is 
highly rarefied to the point that the molecules are so far apart that the 
reflected molecules do not affect the oncoming molecules to any appreciable 
degree. 
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4.4.1 Impingement Pressure: Continuum- Flow Regime 


4. 4. 1,1 The Question 

What impingement pressure formulations are required in continuum- 
flow regime? For an "elemental flat plate" (see Section 4,2.2), 
the most commonly used formulations of impingement pressure (normal 
force per unit area) in this regime are newtonian impact, modified 
newtonian and oblique shock/tangent wedge. 


Newtonian Impact: For newtonian flow, it is assumed that the 

normal component of momentum- is transferred to the plate, while the 
tangential component remains unchanged. The pressure coefficient 
is expressed simply as 

C p = 2 sin 2 a (4. 4. 1-1) 

where a is the impingement angle shown in Figure 4. 4. 1-1. 


Modified Newtonian: Various empirical modifications to Equation 
4. 4. 1-1 have been suggested, to account for the higher pressures 
noted experimentally at small local impingement angles. The 
pressure coefficient is expressed as 

C p = K sin 2 a (4.4. 1-2) 

using a curve fit; e.g.. 


K-C [O.SU + SjM] 

p o a 

where a is expressed- in degrees, and C p 
pressure coefficient". 


o 


is the "stagnation 


(4. 4. 1-3) 


At least three variations of the modified newtonian formulation 
are in use, differing in their treatment of the stagnation pressure 
coefficient: 
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FIGURE 4.4. 1-1 LOCAL IMPINGEMENT ANGLE AT SURFACE ELEMENT 

AND CHANGE IN FLOW DIRECTION FOR NEWTONIAN FLOW 



FIGURE 4. 4. 1-2 IMPINGEMENT CORRECTION FACTOR 
VS LOCAL IMPINGEMENT ANGLE 
IN MODIFIED NEWTONIAN FORMULATION 
AS USED IN "PUMP" (ASSUMING C„ - 1.83) 

r o 
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• Assume a constant (newtonian) value, C = 2. 

p o 

• Assume a constant, empirically-determined value; e.g.. 


• Compute C 



Po 
= P 


from the definition 
o P°° 





1.83. 


(4. 4. 1-4) 


where p Q Is the stagnation pressure (behind a normal shock), 
p OT is the free-stream static pressure, and is the free-, 
stream dynamic pressure. 


Variation of the modified newtonian pressure coefficient with 
impingement angle is shown in Figure 4.4.1 -2, assuming - 1.83. 
This curve was taken from PUMP (Section 3.2.4) 0 

documentation. 


Oblique Shpck/Tangent Wedge: In this method, the impingement 

pressure is assumed to be the same as that existing behind an 
attached oblique shock through which the oncoming flow turns 
parallel to the elemental area. The pressure coefficient is 


C = 3 

P M^-y+i ) 


rfsinVl] 


where 


tan a = 2cot8 


M 2 sin 2 £-1 

M 2 (y+cos2S)+2 


is solved by iteration. A sketch of the oblique shock wave with 
various parameters is shown in Figure 4. 4. 1-3. 


(4.4. 1-5) 


(4. 4. 1-6) 
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FIGURE 4. 4. 1-3 PARAMETERS IN THE OBLIQUE 

SHOCK IMPINGEMENT FORMULATION 
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4. 4. 1.2 Tradeoff Factors/Consi derations 


Implementation: Implementation of- the newtonian impact formulation 

is obvious and straightforward. Implementation of the modified 

newtonian formulation is also straightforward, if a constant Cp is 

assumed. To compute C from its definition (Equation 0 

p o 

4.4.1 -4) will require computation of p Q as a function of y and 
using a small table lookup and standard equations which account for 
the entropy change across a normal shock. This is the approach used 
in PUMP. Implementation of the oblique shock formulation requires 
tabular data for y and throughout the flowfield, as well as an 
iteration scheme to solve Equation 4.4.1 -6. 


"Cost" Factors: It is apparent that the newtonian impact theory gives 

the simplest pressure coefficient form, while the oblique shock 
method is the most complex. We would expect the computer storage 
and execution time requirements to be proportional to the complexity 
of each formulation; iterations may prove very time-consuming. 

Performance Factors : Although the newtonian impact pressure is the 

simplest to implement, its accuracy is typically the worst of the 
three methods considered. The oblique shock and modified newtonian 
methods are comparable in accuracy; within their common domain; 
see Figure 4. 4. 1-4 for representative flat-plate data. However, the 
oblique shock method is applicable only to impingement angles less than 
(approximately) 45°; the exact limiting angle depends upon the local 
flowfield properties. Beyond this limiting angle, the weak, attached 
oblique shock becomes a strong detached shock. Correspondingly, Equation 
4. 4. 1-6 has no solution, so one of the newtonian formulations must be 
used. (The switchover to an alternate formulation must, of course, be 
automatically performed on a per-element basis, as it is in PLIMP.) 


The only method which has the advantage of newtonian impact simplicity 
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PRESSURE COEFFICIENT 



IMPINGEMENT ANGLE, a (DEGREES) 


FIGURE 4.4.1 -4 COMPARATIVE DATA FOR DIFFERENT IMPINGEMENT 
FORMULATIONS (FLAT PLATE WITH ATTACHED 
SHOCK: M ot = 15, T= 1,4) 
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and an expected accuracy comparable to the oblique shock method is 
the modified newtonian formulation. Therefore this formulation is 
considered most suitable for the plume impingement model. 

4. 4. 1.3 Analysis Approach 

To evaluate the' performance of various plume impingement pressure 
formulations, a number of runs of PUMP were made for both MMS/SMM 
and LDEF payloads (Section 3.2.2), for various orientations at an 
on-centerline distance of 100 ft. For all data runs, a mesh size 
of 30 x 30 (Section 4.2.2) was used and object shadowing (Section 
4.2.3) was modelled. 

4. 4. 1.4 Analysis Results 

Table 4. 4. 1-1 provides case-by-case cost and performance data. Note 
that significantly more computer time (roughly a factor of 4) is 
required for the oblique shock method. The computer time required 
by both newtonian formulations is comparable and much less than 
that of the oblique shock formulation. 

In general, the oblique shock and the modified newtonian formulations 
are comparable in accuracy. An exception is case L2, where both 
newtonian methods deviate greatly from the oblique shock result in 
the moment calculation, although the forces are accurate. This is 
because at this orientation, most of the force is contributed by 
the circular end plate, while most of the moment comes from the 
cylinder body. Since the impingement angle on the end plate is 
roughly 70°, the oblique shock option defaults to modified newtonian. 
On the cylinder body, the impingement angle ranges from 20° to grazing 
incidence; the greatest deviations occur for those elements 
immediately neighboring the shadow line. 
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TABLE 4. 4. 1-1 COST AND PERFORMANCE VARIABLES FOR VARIOUS IMPINGEMENT 
PRESSURE FORMULATIONS IN THE CONTINUUM FLOW REGIME 


— 

TEST CASE 

IMPINGEMENT 

PRESSURE* 5 

COST DATA - 

RAW PERFORMANCE DATA 

SENSITIVITY*" 

Number 3 

d 

(ft) 

e 

(deg) 

Execution 

Time 

(Sec) 

Time 

Ratio 0 

F x 

(lb) 

F y 

.(lb) 

H z 

(lb-ft) 

Force Errors 

Moment 

Error 

(+*) 

A Maq 
<+*) 

A Dir 
{+ deg) 

Ml 

100 

20 

OS 

277.0 

- 

.16.71 

5.52 

-17.08 

- 

- 

- 




MN 

69.5 

0.25 

16.77 


-16.81 

0.04 

0.56 

1.58 

M2 

100 

70 

OS 

314.6 

- 

5.35 

0.14 

-14.19 

- 

- 

- 




. MN 

69.6- 

•0.22 

• 5.50 

0.54 

-16.34 

3.23 

4.04 

16.2 

LI 

100 

0 

OS 

131.3 

- 

21.37 

0.0 

0.0 

- 

- 

- 




MN 

48.7 

0.37 

21.37 

0.0 


0.00 

0.00 

- 




NI 

46.8 

0.36 

21.76 



1.82 

0.00 

- 

L2 

100 

20 

OS 

206.9 

- 

18.47 

3.76 

4.18 

- 


- 




MN 

56.4 

0.27 

18.81 

2.84 

6.08 

0.93 

2.92 

45.5 




NI 

42.8 

0.21 

17.74 

3.94 

2.80 

3.59 

1.02 


L3 

100 

70 

OS 1 

252.4 

- 

27.53 

mmm 

29.67 

- 

- 

- 




MN 

40.3 

msm 

27.71 

-8.80 

29.44 

0.24 

0.73 

0.78 




NI 

26.5 

0.11 

25.52 

-8.54 

30.08 

7.22 

0.16 

1.38 

L4 

100 

90 

OS 

115.1 

- 

32.34 

0.0 


- 

- 

- 




MN 

33.4 

0.29 

32.47 



0.40 

0.0 

- 




NI 

27.2 

0.24 

30.86 

0.0 


4.58 

■HHi 

- 


NOTES: 


a M denotes MMS/SMM; L denotes LDEF 

b 0S denotes Oblique Shock 
MN denotes Modified Newtonian 
NI denotes Newtonian Impact 

c With respect to the Oblique Shock result 
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4. 4. 1.5 Preliminary Conclusions and Recommendations 
Based upon the above analysis and results, the modified newtonian 
formulation exhibits the anticipated best compromise between cost 
and accuracy. Therefore this formulation is recoimiended for the 
plume impingement model in the continuum-flow regime. 
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4.4.2 Impingement Pressure; Free Molecule Flow Regime 


4. 4. 2.1 The Question 

What impingement pressure formulations are required in the free-molecule- 
flow regime? For an "elemental flat plate" (see Section 4.2.2), the for- 
mulations of impingement pressure which are usually recommended for the 
free molecule flow regime are newtonian impact, diffuse reflection, and 
thermal accommodation. 


Newtonian Impact : For newtonian flow, it is assumed that the normal com- 
ponent of momentum is transferred to the plate, while the tangential com- 
ponent remains unchanged. The normal and tangential pressure coefficients 


are simply: 

- 2 s o*?© 


(4. 4. 2-1) 


Cp r ® o 


(4. 4. 2-2) 


where 0 is the impingement angle shown in Figure 4. 4. 2-1 (a). 


Diffuse Reflection : Given a stream of incident gas particles, all moving 

with a constant velocity, diffuse reflection occurs when each particle 
hits the plate and is reflected with unchanged energy in a random direc- 
tion; see Figure 4.4. 2-1 (b). The coefficients of pressure are calculated 
in Reference 22 as 


C-pH - 'Z.sCn.Q |su\0 + IT "1 


(4. 4. 2-3) 




(4. 4. 2-4) 


Thermal Accommodation : A rarified gas travels with a stream velocity, c\ 

upon which is superimposed a Maxwellian velocity distribution (Reference 
23) determined from the temperature of the incident gas, T.. . Since the 
motions of the incident gas and the gas reflected from the surface are 
assumed independent (i.e., no collisions), the momentum and energy of 
the reflected gas will not be affected by the presence of the incident 
particles (see Figure 4.4. 2-1 (c)). 
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(a) 


NEWTONIAN IMPACT 




(b). DIFFUSE REFLECTION 



(c) THERMAL ACCOMMODATION 


FIGURE 4. 4.2-1 FREE MOLECULE FLOW IMPINGEMENT FORMULATIONS 
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At the payload surface, incident oarticles undergo collisions with the 
wall molecules. With each collision the particles are deflected by, -and 
exchange energy with, the wall. The momentum and energy of the particles 
leaving the surface can then be defined in terms of two momentum co- 
efficients (a s , a^) and one energy coefficient (a): 


specular accommodation coefficient: ^ ^ 


diffuse accommodation coefficient: 
thermal accommodation coefficient: 


** * 

-Q-TV 


where 


Gt-ew 


X. , P. = tangential and normal momenta of the incident molecules 
T = tangential and normal momenta of the reflected molecules 


w 5 w 


P = tangential and normal momenta of the wall molecules 


E-, E r >' E w = energy of the incident, reflected, and wall molecules. 


At the wall, only a few collisions of each incident molecule while trapped 

within the interstices of the wall are necessary to produce completely 

random emission. However, a much larger number of such encounters is 

necessary to adjust the energy of the incident molecules so that they 

emerge at a temperature consistent with the thermal condition of the 

wall. Experimental data (Reference 11) suggests 0 ^ = = 1 and 0.8 < 

a < 1.0. In subsequent calculations, we assume a = a„ = Ij i.e., the 

d s 

reflected gas leaves randomly from the wall. 

Given the molecule mass (m) , the stream speed (c); and Boltzmann's con- 
stant (k), one can find the molecular speed ratios ($ and S v ): 

s= */(sgi) K 

5 — S' scitO 

The reflected kinetic temperature (T } is then found to be: 

T r = f /£ er< + )]<->-<*> 


(4. 4.2-5) 


(4.4. 2-6) 
(4.4. 2-7) 


(4. 4. 2-8) 
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and the pressure coefficients are: 


** * s *-’ e feGr 1 1 " k*%A 


(4. 4. 2-9) 


CpY — SUr\_©C©S & / \-V -«aJ I 

l XttS v * J 


(4.4.2-10) 


where "erf S v " is the error function of S y . 

When the molecular speed ratio, S, is large, as it normally is in the 
free-molecule region of a rocket exhaust plume, the preceding equations 
can be greatly simplified for most impingement angles. In our situa- 
tion, S will always be greater than ten; so, for angles greater than 
twelve degrees, S y will be greater than two. It is then a good approxi- 
mation to let the error function equal one (erf 2 = .995) and the exponen- 
tial equal zero (e -4 = .02). With these approximations, the simplified 
equations are: 

T r , T- (*2+ S 1 -} (4.4.2-8a) 

*2*0*50 0 + is* + ) (4.4.2-9a) 

“ 2.sor\0co*O (4.4.2-10a) 


Even further simplifications may be reasonable under certain conditions; 
in fact, for S y large enough. Equation 4.4.2-9a reduces to simple newtonian 
impact (Equation 4.4.2-1). 

4. 4. 2. 2 Tradeoff Factors/Considerations 

Imp! ementati on : Implementations of the newtonian impact and diffuse re- 

flection formulations are obvious and straightforward. Implementation 
of the full thermal accommodation formulation requires a subroutine to 
evaluate the error function, and tests for singularities in evaluating 
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the pressure coefficients. 


"Cost" Factors : Newtonian impact theory and diffuse reflection analysis- 

provide the simplest pressure coefficient forms, while the thermal accom- 
modation technique requires somewhat larger computer storage and longer 
execution times. Using the thermal accommodation technique may require 
a table of incident gas temperatures in various regions of the plume 
(depending upon the temperature sensitivity of the pressure coefficients). 

Performance Factors : The newtonian impact pressure coefficients are sim- 

ple in form, but the normal force coefficient is not expected to be very 
accurate, and the tangential force coefficient is zero. 

Classical elastic collisions (treating the gas molecules as hard spheres) 
causing specular and diffuse reflection are often suggested as a mechanism 
for imparting forces to a surface.. However, the conditions for these 
elastic interactions are not satisfied for a fast moving plume gas and 
a typical payload surface (Reference 23). 

The thermal accommodation pressure coefficients are complex but include 
much of the physics of the gas. Further, the thermal accommodation 
technique allows for a "pseudo-diffuse" case when the accommodation 
coefficient (a) is zero. In this case, the energy of the gas is constant 
throughout the interaction and particles leave in random directions from 
the surface. In the degenerate case where the molecular speed ratio (S) 
becomes infinite, the particles travelling to the surface all have a 
uniform speed. Unlike classical diffuse reflection, though, the reflected 
particles are distributed according to a Maxwell -Boltzmann distribution. 

At the other extreme, a = 1, the particles leave the wall with a velocity 

distribution characterized by the wall temperature. Again, experimental 
data suggests 0.8 < a < 1.0. 

The thermal accommodation technique is expected to give the best per- 
formance, in terms of accuracy and generality. 
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4. 4. 2.3 Analysis Approach 


The various free-molecule flow impingement pressure formulations were 
coded on the HPr9825 desk calculator, and a number of plot runs were 
made to compare the pressure coefficients as functions of impingement 
angle. 

Sensitivities of these coefficients to the input parameters were also 
analyzed. 

4. 4. 2. 4 Analysis Results 

Figure 4. 4. 2-2 shows- the normal pressure coefficient as a function of 
impingement angle for various accommodation coefficient inputs using the 
thermal accommodation formulation. At almost normal incidence, the pres- 
sure coefficient is rather sensitive to the value of the accommodation 
coefficient, varying by almost a factor of 2 over the input accommodation 
coefficient range of 0 to 1. 

Figure 4. 4. 2-3 shows the variation of the normal pressure coefficient as 
a function of impingement angle for various incident gas temperatures. 

The pressure coefficient is not very sensitive to the incident gas tempera- 
ture even at very low impingement angles, where a factor of 4 increase 
in temperature only increases the normal pressure coefficient by 10%. 

Figure 4.4. 2-4 depicts the variation of the normal pressure coefficient 
as a function of impingement angle for various wall temperatures. Over 
the range of pay.load surface temperatures expected in space, the pressure 
coefficient varies by 2-5%. 

Figure 4.4. 2-5 shows the variation of the tangential pressure coefficient 
as a function of impingement angle for various incident gas temperatures 
(note from Equation 4.4.2-10 that the accommodation coefficient and wall 
temperature do not affect the value of the tangential pressure coefficient). 
Only at small impingement angles does the pressure coefficient differ from 
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FIGURE 4. 4. 2-5 VARIATION OF Cp T /SIN-0 COS 0 AS A FUNCTION 
OF 0 FOR VARIOUS INCIDENT GAS TEMPERATURES 
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the value of 2,0 predicted by diffuse reflection. (In newtonian flow, 
the tangential pressure coefficient is always zero.) 

A comparison of various impingement recipes considered for free-molecule 
flow is given in Figure 4. 4. 2-6. Here, the normal pressure coefficient 
is plotted as a function of impingement angle using the newtonian impact, 
diffuse reflection, and thermal accommodation formulations (a = 0, 0.8, 
0.9, 1.0), Also plotted, for comparison only, is the recommended impinge- 
ment formulation in the continuum regime: modified newtonian impact (see 

Section 4.4.1 ) . 


A comparison of storage and execution times among the free-molecule impinge- 
ment formulations using the HP-9825 showed that the thermal accommodation 
formulation took an average of 95 msec to execute and about 500 storage 
registers. The other simpler formulations were executed in about 18 msec 
using 80 storage registers. 


From plume information documented in Reference 24, the value of the mole- 
cular speed ratio (S) in Equation 4. 4. 2-6 is about 11.5. For all 
impingement angles larger than seven degrees, an error of less than 1% is 

-S 2 

incurred by using the approximations erf S y = 1 and e v = 0 (Equations 
4.4.2-8a, 9a, and 10a). These approximations may speed up execution time 
by a factor of two. 


4. 4. 2. 5 Preliminary Conclusions and Recommendations 

There is a large variation in the values of the normal and tangential pres- 
sure coefficients in the free molecule regime, depending on which formula- 
tion is used; further, none of these formulations agrees very closely with 
the modified newtonian impact formulation used in the continuum regime. 

From our present understanding of the physics of the gas/surface inter- 
action, we believe the thermal accommodation formulation to be the 'most 
accurate. 


Newtonian impact, theory has been seen to give at best a 10% underpredic- 
tion in the normal pressure coefficient, and does not predict any tangen- 
tial forces. The diffuse reflection formulation overpredicts the normal 
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FIGURE 4.4. 2-6 COMPARISON OF C pH FOR VARIOUS IMPINGEMENT FORMULATIONS 
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forces by 30%, but accurately predicts tangential forces, except at shallow 
impingement angles. 

The thermal accommodation formulation generates pressure coefficients 
which are very sensitive to the value used for the accommodation- co- 
efficient, with conservative force calculations using smaller values of 
the accommodation coefficient. These coefficients are fairly insensitive 
to input values of gas temperature and wall temperature. To within 1 % 
error, the error function can be set to unity and the exponential to zero 
in these pressure coefficients for impingement angles greater than seven 
degrees, when considering RCS plumes. 

* 

At small impingement angles, tangential pressure coefficients are not 
small; in fact, at angles as large as 30 degrees, tangential forces are 
larger than normal forces. 

We therefore recommend: 

• Use of the thermal accommodation pressure coefficient formulation 
« Use of a constant incident gas temperature 

o 

• Use of the approximations erf-S = 1 and exp (-$ v )‘= 0 for 
impingement angles greater than seven degrees. 
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4.4.3 Impingement Pressure: Transitional -Flow Regime 


4.4. 3.1 The Question 

What impingement pressure formulation is required in the transitional 
flow regime? This regime is defined in terms of the local knudsen 
number, Kn: 0.01 < Kn < 10 for transitional flow. In this regime, 
direct computation of impingement pressure by analysis of flow 
mechanisms is extremely complex (see Section 4.4). Therefore, the 
established practice is to use a "weighting function" to transfer 
from continuum to free molecule impingement formulations. This weight- 
ing function, f(Kn), varies from 0 to 1 as the knudsen number varies 
from .01 to 10. Then the normal and tangential pressure coefficients 
in the transitional regime are computed as: 

C V C '’w'* 0c " j v 

(4 

where "tf" denotes transitional flow 

"fmf" denotes free-molecule flow, 

"cf" denotes continuum flow 

Both "hard" and "soft" transitions between flow regimes will be con- 
sidered in this study. 


Hard Transition : The simplest weighting function is a step function. 
Letting the step transition occur at the "logarithmic midpoint" of 
the knudsen number interval , 



V* <0- 
VCw 


(4.4. 3-2) 


Soft Transitions : These transitions use weighting functions which are 
continuous over the knudsen number interval. As a first approxima- 
tion, let f(Kn) vary linearly with Log Kn. Then, as can be seen 
in Figure 4.4. 3-1, 




(4.4. 3-3) 
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FIGURE 4.4.3-1 A FIRST-APPROXIMATION SOFT-TRANSITION WEIGHTING FUNCTION 
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The PLIMP program (Section 3.2.4) uses a smooth "s-curve" weighting 
function, which also has a continuous first derivative. This weight- 
ing function varies as the square of a sine function: 


few* *** 

Finally., in order to reduce execution time. Equation 4. 4. 3-3 can 
be replaced by a more accurate linear weighting function based on 
a simple expansion (slope match) of Equation 4. 4. 3-4: 

In effect, the transition flow limits, 10"^ 5 Kn £ 10, are changed 
to the function breakpoints, .035 < Kn <2.85. The weighting func- 
tions listed above are shown in Figure 4. 4. 3-2. 

4.4. 3.2 Tradeoff Factors/Considerations 

Implementati on : Implementation of a transition regime impingement 

formulation will require the ability to compute both continuum and 
free-molecule flow pressure coefficients, a routine and appropriate 
flowfield data to calculate the local mean-free-path , and specifica- 
tion of a characteristic payload subshape dimension to calculate the 
knudsen number. With these data, any of the previous weighting func- 
tions will give pressure coefficients for the transitional regime. 

“Cost 11 Factors : The knudsen number calculation has the largest impact 
on storage requirements, but a minor impact on execution time, be- 
cause it is only done once for each subshape, while the impingement 
pressure is computed for each element (see Section 4.2.2). 

The hard transition obviously requires the minimum storage and ex- 
ecution time, since the weighting function computation reduces to a 


(4.4. 3-4) 


(4.4. 3-5) 
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FIGURE 4. 4. 3-2 WEIGHTING FUNCTIONS CONSIDERED IN THE STUDY 








single test, and only one set of pressure coefficients is computed. 

The slope-match linear fit of the sin function is the fastest of 

the soft transitions, since it "shrinks" the region for which dual 

2 

pressure coefficients need be computed. The sin transition is the 
only case where the time for computation of the weighting function 
itself is comparable to the time to compute the two sets of pressure 
coefficients. 

o 

Performance Factors : The soft sin transition is considered the 
most accurate form, and has a continuous first derivative. The hard 
transition obviously has the largest error, which occurs in the 
middle of the transitional region.- However, if the continuum and free 
molecule pressure coefficients are close in value, then the transi- 
tional coefficients will be relatively insensitive to the weighting 
function. The slope-match linear fit is not expected to vary much 

9 

from the sin form, and may be preferred when execution time is 
criti cal. 

4. 4. 3. 3 Analysis Approach 

The performance of the hard transition and the soft sin 2 weighting 
function, used with Equation 4. 4. 3-1, were evaluated by a series of 
plot runs on a desk calculator. The reconmended continuum and free 
molecule pressure coefficients (modified newtonian and thermal ac- 
comodation, respectively; see Sections 4.4.1 and 4.4.2), were used 
in Equation 4.4. 3-1 and transitional pressure coefficients were 
calculated as a function of the knudsen mmber for various impinge- 
ment angles. A value of 0.9 was used for the acconmodation coefficient. 

4.4. 3-4 Analysis Results 

Figures 4. 4. 3-3 through 4. 4. 3-7 show the transition from continuum 
pressure coefficients (normal and tangential) to free molecule 

O 

pressure coefficients using hard and soft (sin ) weighting functions. 
Assuming the sin transition to be the most realistic, the maximum 
error incurred by using the hard transition occurs at Kn B 0.316 
as expected . 


no 





FIGURE 4. 4. 3-3 HARD VS. SOFT TRANSITION (SIN 2 FORM): 0 = 90° (NORMAL INCIDENCE) 




FIGURE 4.4. 3-4 HARD VS. SOFT TRANSITION <SIN Z FORM): 0 - 60° 
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FIGURE 4.4. 3-5 HARD VS 
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FIGURE 4.4. 3-6 HARO VS 


tangential pressure coefficient, C, 
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FIGURE 4.4. 3-7 HARD VS. SOFT 


tangential pressure coefficient, 



knudsen number, Kn 


(SIN 2 FORM): 6 - 0° (GRAZING INCIDENCE) 



For the normal pressure coefficient, the maximum error is about 15%. 

The error is fairly insensitive to impingement angle, but does in- 
crease somewhat at shallow impingement angles. Since the tangential 
pressure coefficient in the continuum regime is always zero, the hard- 
transition value of Cpj is exactly twice the soft-transi tion value 
• at Kn = 0.316, and the maximum error is always 100%. 

O * 

To consider the linear expansion of the soft sin transition (Eqn 4,4. 3-5), 
the midpoint tangent to the soft transition curves in Figures 4. 4. 3-3 
through 4. 4. 3-7 can be drawn. The error is largest at the breakpoints, 

Kn = .035 and Kn = 2.85, and is typically on the order of 2 or 3 percent. 

4.4. 3-5 Preliminary Conclusions and Recommendations 

Based on the above analysis, the hard transition formulation is not 

recommended for determining impingement pressure coefficients in 

the transitional regime. The soft sin formulation is simple and 

o 

accurate, but the linear expansion of the sin formulation is quicker 
with comparable accuracy. Therefore, this linear weighting func- 
tion is recoranended. 
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4.5 Dynamical Effects of Modelling Alternatives 

In an on-orbit simulation incorporating a plume impingement model, 

each thruster firing induces forces and moments on an object (payload or 

RMS arm) immersed in the plume, thus exciting the dynamics of the object 

and perturbing its state. In general, these objects are also subjected 

to non-impingement forces and moments, which also perturb the object 

state. 

The presence of substantial non-impingement disturbances tends to reduce 
the importance of plume-modelling accuracy for accurate determination of the 
object state (which is, of course, the ultimate goal of the simulation). 

The nature of the object dynamics may either "amplify" or "attenuate" state 
errors due to plume-modelling inaccuracies. 

The studies documented in the following subsections were conducted to 
explore relationships between plume-impingement disturbances, non- 
impingement disturbances, and object dynamics, to assess the overall 
importance of plume modelling in computing an accurate object state. 
Specifically: 

$ Section 4.5.1 considers plume impingement on the RMS, a situation for 
which impingement disturbances may be small relative to other 
disturbances. 

• Section 4.5.2 considers propagation of plume disturbances through 
the dynamics of a free-flying payload, to determine whether 
plume-modelling errors are amplified by the vehicle dynamics. 
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4.5.1 RMS Arm: Impingement Moments vs. Inertial -Reaction Moments 

4. 5. 1.1 The Question 

Are plume impingement forces and moments on an RMS arm significant in 
comparison to other forces and moments acting on the RMS? One non- 
impingement effect acting on the RMS, even when empty, is the "inertial 
reaction" experienced whenever a thruster firing imparts linear and 
angular accelerations to the Orbiter. 

If a thruster fires while the RMS bears an attached payload, the RMS will 
experience even greater inertial -reaction effects. In addition, there 
may be plume impingement on the payload, which will also induce forces 
and moments on the RMS. 

If it becomes obvious that impingement effects are insignificant compared 
to these other effects, then impingement calculations for the RMS itself 
need not be incorporated in the new plume model. 

As a "by-product" of this sensitivity analysis, we hope to determine 
whether either impingement or inertial moments (or both) will approach 
RMS design limits for the maximum moment allowed at the shoulder. 

4. 5. 1.2 Tradeoff Factors /Considerations 

Implementati on : In order to provide for RMS plume impingement analysis 

in the plume model, the RMS simulation must provide the position and 
orientation of each segment of the arm. The plume model must then cal- 
culate forces on each segment and moments about each joint. The RMS 
simulation must then accept these force/moment data for use in dynamical 
calculations on actuator responses, etc. 

Cost : The cost will be large. The plume model itself will require RMS 
geometry parameters as well as "multiple object" storage and logic, to 
be able to calculate impingement effects on a payload and an RMS arm 
simultaneously. An impingement force on a single RMS segment will create 
moments about all RMS joints. Execution time will increase substantially 
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as impingement is determined for each segment of the RMS. A large in- 
crease in execution time is also expected in the RMS dynamics model, 
which must accept the force/moment data and compute updated RMS positions. 

Performance ; Incorporation of RCS impingement calculations on the RMS 
will require greater accuracy in the continuum flowfield regime, and 
precise modelling of Orbiter surface shadowing. The improvement in 
overall RMS modelling accuracy is expected to be small, for two reasons: 

• With an attached payload, the forces and moments due to inertial re- 
actions and impingement on the payload will be vastly greater than 
those due to impingement on the RMS itself,. 

• Even for an empty arm, RMS dynamics, being nonlinear, will be affected 
perceptibly only when net moments are comparable to the “saturation" 
limits of RMS mechanical components - actuators, brakes, and clutches. 

4.5. 1.3 Analysis Approach 

Theory : The relationship between “fictitious forces" (e.g., centrifugal 

and coriolis forces) as seen in a non-inertial reference frame, and 
forces seen in a fixed inertial frame, is found in many elementary 
physics books. Using the position vectors in Figure 4. 5. 1-1, the torque 
about the RMS shoulder as measured in the non-inertial Orbiter frame, 
due to thruster forces, is: 


r s * (T. rn L tWj Z t )x(2T x^) - u?(T Si.*) + t *>.Z £ (t. •T)(4.5.1-l) 

where parameters are as defined in the figure. 


Computati ons : In this analysis, only the forward +Z RCS thruster was 

considered. Impingement effects on the RMS are expected to be largest 
for this thruster. Plume impingement moments were estimated by using 
PLIMP (Section 3.2.4), including Orbiter shadowing, and using a modified 
newtonian impingement formulation. 
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Parameters 

F a, force 

r, r , Z ^ position vectors 


m a, mass 

I a, moment of inertia 
a a, angular acceleration 


of rotating frame 


T Aj moment 
Subscri pts 

o a, Orbiter 
f a. Fixed (Inertial) 

■j ^ i th element of RMS 



^ if. * F. 

OJ - Z 


FIGURE 4 . 5 . 1-1 


LES INVOLVED IN CALCULATION’ 
S ABOUT RMS SHOULDER 


OF INERTIAL-REACTION 
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Inertial and impingement moments were determined for two basic orienta- 
tions of the RMS. For the first orientation, the arm was considered to 
be a straight uniform rod with' a mass of about 26 slugs (830 1b.), fifty 
feet in length, at an angle 0 to the longeron, as depicted in Figure 
4.5. 1-2. This case (with 9 = 90°) was chosen since it is expected to 
have large inertial and impingement forces. As the arm is folded, im- 
pingement forces are expected to' decrease faster than inertial forces, 
although a bent arm is operationally more typical. For the second orien- 
tation (the "maximum impingement" orientation), the arm is bent forward 
over the cabin, and passes above the center forward +Z RCS thruster; see 
Figure 4. 5. 1-3. 

Equation 4. 5. 1-1 assumes that the initial angular velocity was zero; the 
Orbiter was considered to be initially inertially stabilized. The RMS 
was assumed to be rigidly attached to the Orbiter, and only the moment 
about the shoulder was computed. 

4.5. 1.4 Analysis Results 


The inertial moment for the straight arm case was found to be 



- 3o6eos0 +3S%(M.Q 



(4. 5. 1-2) 


For the case depicted in Figure 4.5. 1-2, where 9 = 90°, the moment is 

235 lb-ft, compared to an impingement moment of about 55 lb-ft found with PUMP. 

The three dimensional moment results are given in Table 4.5. 1-1. 


Moments found in the maximum impingement orientation exceed the RMS 
specification of a guaranteed 772 lb-ft stall torque. 

The RMS will slip, and may possibly be damaged. This orientation is 
not considered typical of normal operations, but shows that impingement 
forces may place operational constraints on the RMS, 
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FIGURE 4.5.1-; 



STRAIGHT RMS ORIENTATION 
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4. 5. 1.5 Preliminary Conclusions and Recommendations 


Based upon the above analysis-, we have reached the following conclusions 

• Simulation of plume impingement on the RMS will make little or no 
contribution to the accuracy of RMS dynamical analyses, for normal 
operations. 

• Plume impingement considerations can be used to define gross opera- 
tional constraints on the RMS (such as keeping the arm away from the 
vicinity of RCS thrusters). 

We therefore make the following recommendations: 

• Impingement calculations on the RMS arm should not be incorporated 
in the plume model used in integrated simulations. 

• Offline "stand-alone" analysis programs should be used to define 
plume-related constraints on the RMS. 
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4.5.2 Free Flyers: Propagation Effects 


4. 5. 2.1 The Question 

Do propagation effects "amplify" plume modelling errors? Errors in geo- 
metric representation, flowfields, pressure coefficients, etc. all con- 
tribute to errors in the force and moment computations. These errors in 
turn cause errors in the angular and linear motion computed in an inte- 
grated simulation. Although these errors may initially be small, they 
tend to grow with time for two reasons: 

• Initial errors are propagated through the vehicle dynamics. 

• At any later time, forces and moments will be computed using a state 
which is in error (relative to the true state; see Figure 4. 5. 2-1), 
thus "compounding" the error. 

There was some concern that under some circumstances, these effects might 
cause state errors much larger, percentagewise, than the underlying, plume 
modelling error— say, a 20% error in the plume model causing a factor^of- 
two error in the state perturbations. 

i 

4. 5. 2. 2 Tradeoff Factors /Considerations 


If we find significant amplification in the errors, the plume modelling 
errors will have to be kept low to have reliable results in the integrated 
simulation. This will require higher fidelity and more complex modelling. 

Implementation : To reduce plume modelling errors, the flowfield and im- 
pingement recipes must be computed from complex equations and logic. Finer 
mesh sizes and more accurate payload geometry representation must be in- 
cluded in the program. 

Cost : The cost will be -high in terms of complexity, storage, and execu- 

tion time if plume modelling errors are to be held to a very low level. 

Performance : By knowing whether propagation effects amplify plime model- 

ling errors, accuracy requirements for the new model can be derived. Long 
time duration approaches, and approaches allowing intense plume impinge- 
ment, may require more accurate plume modelling, since the propagation 
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FIGURE 4. 5. 2-1 PROPAGATION ' EFFECTS ON PLUME MODELLING ERRORS. 
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effect will become more visible in these cases. 


4. 5. 2. 3 



Theory : In this sensitivity study, we will only be concerned with plume 

modelling errors which scale the force on an elemental area by a con- 

■+ 

stant, a. The elemental force dF is formed from the elemental area dS, • 
the dynamic pressure q, and the pressure coefficient Cpi 

dF ^ Of ^ ^ (4. 5. 2-1) 

For a plume model where the dynamic pressure or the pressure coefficient 
has such a constant error, then 


dT - oi C p dS 

The force and moment calculations on a payload become scaled by a, causing 
instantaneous translational and rotational acceleration errors which de- 
pend linearly on a. These acceleration errors in turn cause payload state 
vector errors. 

Computation : The Quick Look Disturbance Program, QLDP, (Section 3.2.4) 

was modified to determine propagation effects due to a constant plume 
modelling error. Linearized equations of motion (sometimes called the 
Clohessy-Wiltshire equations) were added for realistic propagation of the 
payload state. Various approaches to the gravity-gradient stabilized 
LDEF (Section 3.2.2) were selected from Shuttle Engineering Simulator 
(SES) data: a nominal braking approach (Reference 25), a V-bar approach 

using X and Z jets, a V-bar approach using only X jets, and an R-^bar 
approach (Reference 26). For each approach, the payload state was deter- 
mined using constant scaling factors of 0.50, 0.90, 0.95, 1.00, 

1.05, 1.10, and 2.00. The + 10% variation (.90<a<1.10) is considered 
typical of the accuracy sought in the new model. Factor-of-two variations 
are considered typical of simplified models (e.g. , see Reference 10), 
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4. 5. 2. 4 Analysis Results 


The final payload states (position, pitch, and velocities) at the end of 
the various approaches are given in Tables 4. 5. 2-1 through 4. 5, 2-4 for 
various scale factors. Also shown in these tables are the ratios of the 
payload state variables to the payload states in the nominal case (con- 
stant a=l). These ratios show that in most cases, constant errors in 
the plume modeX translate directly into errors in the final payload state 
with no "amplification" at least for LDEF. 

To indicate the significance of such state errors for man-in-loop simula- 
tion, Figures 4. 5. 2-2 through 4. 5.2-5 show how the LDEF position and 
orientation in the local vertical coordinate system would look at the end 
of each approach. For the approaches using +2 thrusters, the translational 
displacement is noticeably affected by the scale factor a, although the 
pitch is not significantly different. Shuttle approaches with the +Z 
thrusters inhibited show that "a small error in a small perturbation is 
very small". 

4.5. 2.5 Preliminary Conclusions 

Based on the above analysis, constant errors in the force determination in 
a plume model translate directly into errors in our knowledge of the final 
payload state. Propagation effects do not in general amplify plume model- 
ling errors. 
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TABLE 4.5. 2-1 FINAL LDEF STATE VS. SCALE FACTOR 
FOR A NOMINAL BRAKING APPROACH 


SCALE FACTOR 

STATE PERTURBATIONS 

RATIOS OF PERTURBATIONS ; 

a 

AO 

AX 

AZ 

AG 

Av x 

!!■ 

AG 

Ar 

AG ' 

Av 


(deg) 

(ft) 

(ft) 

(deg/sec) 

(ft/sec) 


A0„„ m 

nom. 

Ar nom. 

^nom. 

Av nom. 

0.50 

-4.468 

-16.23 

-6.23 

-.0244 

-.1408 

-.0175. 

.51 

.51 

.50 

.53 

mmmm 

-7.960 

-28.46 

-11.29 

-.0441 


-.0353 

.90 

.90 

.90 

.95 

KB 

-8.388 

-29.95 

-11.93 

-.0465 


-.0377 

.95 

.95 

.95 

.95 

1.00 (nominal) 

-8.815 


-12.56 

-.0489 

-.2642 ' 

-.0401 

1.00 


1.00 

1.00 


-9.239 


Hi 


-.2757 


1.05 

1.05 

1.05 

1.04 

IBB 

-9.661 

v . 

m 


-.2872 ■ 


1.10 

1.09 

1.10 

1.09 


lliggll 

-59.72 

-25.04 

“.0935 

-.4745 

-.0909 

1.92 

1.91 

1.91 

1.81 
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TABLE 4. 5. 2-2- FINAL LDEF STATE VS. SCALE FACTOR FOR A 
Y-BAR APPROACH USING +Z THRUSTERS 


SCALE FACTOR 

STATE PERTURBATIONS 

RATIOS OF PERTURBATIONS 

a 

A0 

AX 

AZ 

A0 

t> 

< 

X 

BB 

A0 

Ar 

A0- ‘ 

Av 


(deg) 

(ft) 

(ft) 

(deg/sec) 

(ft/sec) 



Ar nom. 

A ®nom. 

Av nom. 

0.50 

0.1718 

-18.46 

18.81 

.0060 

-.0181 

Bsa 

.82 

.52 

♦ 

*>4 

.52 

mm 

0.1916 

-31.86 

33.37 

.0076 

-.0275 

.1077 

.91 

.90 

.97 

.91 

Di 

0.1996 

-33.46 

35.17 


-.0285 

.1133 

.95 ' 

.95 

.99 

.95 

1.00 (nominal) 

0.2096 

-35.05 

36.96 


BHi 

.1189 

1.00 

1.00 

1.00 

1.00 

1.05 

H 

-36.63 

38.74 

■ 

■ 

m 

1.06 

1.05 

1.01 

1.05 

1.10 


-38;l9 

40.52 

BIB 



1.13 

1.10 


1.09 

2.00 

0.8920 

-64.29 

71.68 

.0101 

-♦0412 

.2242; 

4.26 

1.89 

1.29 

1.86 
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TABLE 4.5. 2-3 FINAL-LDEF STATE VS. SCALE FACTOR FOR A 

V-BAR APPROACH WITH +Z THRUSTERS INHIBITED 


SCALE FACTOR 

STATE PERTURBATIONS 

RATIOS OF PERTURBATIONS 

a 

A0 

AX 

AZ 

AO 

4v x 

AV Z 

to 

Ar 

A0 

Ay 

! 


(deg) 

(ft) . 

(ft) 

(deg/sed 

(ft/sec) 

(ft/sec) 

A0 " 
nom. 

Ar 

nom. 

A ®nom. 

Av '1 

nom. 

0.50- 

0.00816 ' 

-.046 

1056 

-.00003 

.00001 

.00016 

.54 

.50 

.50 



..50 

0.90 

0.01466 

-.084 

1101 

-.00005 

.00002 

.00028 

.97 

.90 


.88 

0.95 

_ 0.01438 

-.089 

.106 

-.00005 

.00002 

.00030 

.95 

.85 

1 

.94 

1.00 (nominal) 

' 0.01513 

-.094 

.112 

-.00006 

.00002 

.00032 

1.00 

1.00 

1.00 

1.00 

1.05 ' 

0.01588 

-.099 

.118 

-.00006 

.00002 

.00033 

1.05 

' 1.05 

1.00 ' 

1.05 

1.10 

, 0.01663 

-.103 

.123 

-.00006 

.00002 

' .00035 

1.10 

1.10 

© 

1.00 

1.11 

2.00 

0.03014' 

-.188' 

.224 

-.00011 

.00004 

.00063 

1.99 

2.00 

1.83 

1.97 



%i.V 
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•TABLE 4. 5.2-4 FINAL LDEF STATE VS. SCALE FACTOR 
FOR AN R-BAR APPROACH 


SCALE FACTOR 

STATE PERTURBATIONS 

RATIOS OF PERTURBATIONS , 

a 

A0 

AX 

az 

AO 

4v x 

AV Z 

AO 

Ar 

AQ ' 

Av i 


(deg) 

(ft) 

(ft) 

(deg/sed 

(ft/sec) 

(ft/sec) 

A0__ 

non. 

Ar- „ 
non. 

A ®nom. 

Av * 

nom. 

0.50 

.00062 

-.050 ' 

-.025 

<1Q" 5 

-.00005 

-.00001 

.50 

.50 

m 

.46 i 

i 

0.90 

.00111 

-.089 

-.046 

<10" 5 

-.00010 

r.COGOl 

.90 

.89 

■■■ 

.91 

0.95 

.00117 

-.094 

-.048 

<10‘ 5 

-.00010 

-.00001 

.95 

.94 

biiik 

.91 

1.00 (nominal) 

.00123 

-.100 

-.051 

<10- 5 

-.00011 

-.00001 

1.00 

1.00 

1.00 

1.00 

1.05 ' 

. .00129 

-.105 

-.054 

<10“ 5 

-.00011 

-.00001 

1.05 

1.05 


— 

1.10 

.00136 

-.110 

-.056 

<10" 5 

-.00011 

-.00002 

1.10 

1.10 

bSSI 

IB 

2.00 

. .00246 

-.199 

-.102 

<10“ 5 

-.00021 

-.00002 

2.00 

2.00 


1.91 
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FIGURE 4. 5.2-3 FINAL LDEF STATE AS A FUNCTION OF SCALE FACTOR 
FOR A V-BAR APPROACH USING +Z THRUSTERS 
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FIGURE 4. 5. 2-4 FINAL LDEF STATE AS A FUNCTION OF SCALE FACTOR 
FOR A V-BAR APPROACH WITH +Z THRUSTERS INHIBITED 


1 67 




4.6 Contaminant Flux Modelling Approaches 
(This section to be provided at a later date). 
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5. CONCLUSIONS AND RECOMMENDATIONS 


In this section, the preliminary conclusions' from the individual 
sensitivity studies are integrated. These preliminary con- 
clusions are compared, and potential interactions are considered. 
Based on these interactions and intended model uses, a final set 
of conclusions and recommendations is determined. 

5.1 Integration of Individual Study Results 

Results obtained in the individual sensitivity studies of Section 4 
were based in part on assumptions concerning the outcome of other 
individual studies (e.g., the shadowing studies assumed radial 
streamlines). In general, though, "optimum" conditions were estab- 
lished for all factors except the factor under study;, therefore 
few, if any, of the preliminary conclusions were expected to change. 

5.1.1 Orbiter/Jet Geometry Studies 

Analysis of Orbiter shadowing of RCS plumes (Section 4.1.1) was 
done assuming radial flow streamlines and axisymmetric nozzles. 

The results indicated that Orbiter shadowing should be- included 
in the model. Actually, many Orbiter thrusters have some scarfing; 
further, the flow streamlines are not precisely radial in the imme- 
diate vicinity of the thrusters. However, there is no reason to 
believe that either of these effects would change the preliminary 
conclusion of this study. 

The analysis of Orbiter reflection (Section 4.1.2) was based upon 
the flow of an aft +Z engine interacting with the Orbiter tail. 

The conclusion of this study was that no reflection away from 
Orbiter surfaces toward a payload occurs. Although the details of 
flow interaction might be different for other engines and other 
Orbiter surfaces, the basic conclusion' is expected to hold in all 
cases. 
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5 J ;2 Object Geometry Studies 


The study of object geometry representation indicated that a drastic 
simplification of object geometry yields unacceptable accuracy, but 
a "first-order" simplification (e.g., treating a twelve-sided poly- 
gonal prism as a circular cylinder) is desirable. The latter con- 
clusion obviously depends on the element subdivision mesh size used 
in the analysis. The mesh size analysis (Section 4.2.2) showed 
that a 10 x 10 mesh gave almost the same results as a 50 x 50. Hence, 
one would expect a 12 x 12 cylinder (which in fact is analyzed as a 
twelve-sided polygonal prism) to have essentially the same impinge- 
ment as a circular cylinder. Therefore, mesh size considerations do 
not change the preliminary conclusion, but in fact add support. 

Assumptions made in the mesh size analysis are justified since various 
geometries were studied in various regions of the plume using a so- 
phisticated flowfield. 

A recommendation was also made to Include payload shadowing (Section 
4.2.3). This recommendation was based on the results of PL'IMP analyses 
of both LDEF and MMS-SMM, using fine mesh sizes. Again, the pre- 
liminary conclusions are not changed. 

5.1.3 Flowfield Studies 

Conclusions reached in the flowfield analyses are: (.1) both scarfed 
and unscarfed RCS thrusters (Section 4.3.2) can be identically mo- 
delled in terms of a patched source flow and radial streamlines 
(Section 4.3.1); (2) the ambient atmosphere has no apparent effect 
on the plume (Section 4.3.3); and (3) implementation of a "practical 
boundary" test will save computer time without meaningful sacrifice 
of accuracy for operations of interest (Section 4.3.3). These re- 
sults are independent of any assumptions concerning the payload or 
Orbiter; consequently the conclusions remain unchanged. 

One conditional conclusion was also reached in the flowfield studies: 
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an "equivalent engine" formulation should be implemented, if simul- 
taneous firings of adjacent engines occur frequently in operations’ 
of interest (Section 4.3.4). Additional input from NASA will be 
required to arrive at a final decision on implementing this feature. 

5.1.4 Impingement-Pressure Studies 

Analysis of impingement in the continuum regime (Section 4.4.1) sug- 
gested that a modified newtonian formulation be used. This conclu- ■ 
sion is based upon a PLIMP analysis of LDEF and MMS-SMM payloads at 
various orientations. Payload shadowing and a fine mesh size (30 x 30) 
were assumed to be important and were, therefore, included. 

The free-molecule impingement analysis (Section 4.4.2) was performed- 
analytically without assumptions regarding payload, mesh size, etc. 
Results of other studies cannot, therefore, affect the free-molecule 
impingement conclusion— to use a thermal -accommodation formulation. 

The transitional impingement analysis (Section 4.4.3) was also analy- 
tic and based upon impingement results in the continuum and free- 

molecule regimes. Use of a "soft" linear weighting function match- 

2 

mg the midpoint slope of the sin "S-curve". was 1 recommended. 

Since these impingement results are based on valid assumptions, no 
change is recommended in the preliminary conclusions. 

5.1.5 Dynamical Effects Studies 


In considering dynamical effects of modelling alternatives, the pre- 
liminary conclusions in Section 4.5.1 were that impingement on the 
RMS arm should not be modelled; rather, an offline "stand-alone" 
program should be used to define gross operating constraints on the 
RMS. These conclusions were based on PLIMP analysis of the RMS, 
modelled as a set of cylindrical segments, impinged upon by an axi- 
symmetric thruster, with Orbiter shadowing included. Inertial re- 
actions of the RMS to RCS firings were analytically determined. 
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Since all the PLIMP assumptions have been determined to be good, no 
change to the preliminary conclusions is suggested. 


For free-flying payloads (Section 4.5.2), constant ‘'scale factor" 
errors in the plume pressure model translated directly into errors 
in knowledge of the payload position and orientation (i.e., no "am- 
plification" effects), in most cases. This result, however, was 
based upon an analysis where the LDEF payload was modelled as four 
plates (not a first-order simplification) and the impingement was 
simple newtonian. Since these have been seen not to be good assump- 
tions, and a scale factor is not the most general form of modelling 
error, the conclusion of this study must be relaxed. Still, we . 
expect no significant amplification effects from vehicle dynamics. 

5,1.6 Contamination Studies 


(This section will be supplied at a later date.) 

5.2 Summary of Modelling Recommendations 

Based on the results obtained in the Plume Impingement Model Sensi- 
tivity Analysis Study, the recommendations of Table 5-1 are suggested 
for implementation in the Level B requirements for the plume model. 
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TABLE 5-1 SUMMARY OF MODELLING RECOMMENDATIONS 


MODELLING AREA 

RECOMMENDED IMPLEMENTATIONS 

4.1 

Orbiter/Jet Geometry 

• Include shadowing caused by the Orbiter moldline, cabin, OMS-pod bulkhead, and tail. 

• Do not consider reflections of the plume. 

4.2 

Object Geometry 

• Include some "first order" simplification in the object modelling. 

• The mesh size should be variable, and on the order of 10x10 in the near field. 

• Include back-side shadowing, and shadowing of one payload element by another payload subshape. 

4.3 

Flowfield 

• Model the flowfield as patched source flow with radial streamlines. 

• Neglect effects on the flowfield due to nozzle scarfing and ambient atmosphere effects. 

• Include a "practical boundary" test. 

• If simultaneous firings of adjacent engines occur frequently, implement an\"equ1valent engine" 

formulation. ... 

4.4 

Impingement 

• Calculate impingement pressure in the continuum regime using a modified newtonian formulation. 

• Calculate impingement pressure in the free-molecule regime using thermal accommodation analysis; 
use approximations for large impingement angles. 

t In the transitional regime, use a linear weighting function, based on the sin 2 function, to 
combine continuum and free-molecule results. 

4.5 

Dynamical 

Effects 

• 

• Neglect plume impingement on the RMS for integrated simulations. Use off-line "stand-alone" 
analysis to Identify gross RMS operating constraints. 

4.6 

Contamination 

To be supplied. 


Additional NASA Input required for final decision. 
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